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As a consequence of the world energy supply problems large kilo- 
watt and megawatt capacity wind turbine generators are being designed 
and built to take advantage of the energy in the wind. These appara- 
tus are equipped with synchronous alternators for direct synchroniza- 
tion with utility networks. Because a wind turbine generator inher- 
ently has unique power system features - such as a large moment of 
inertia relative to its power capacity and a continuously varying in- 
put torque - its electrical behavior must be investigated ^either by 
testing or by analysis. The transient response of the system to elec- 
trical faults is one mode of the system behavior that must-be studied 
analytically rather than experimentally. 

In order to obtain a measure of its response to short circuits 
a large horizontal axis wind turbine generator is modeled and its per- 
formance is simulated on a digital computer. The baseline model for 
the simulation is the DOE/NASA 100 kW wind turbine generator located 
near Sandusky, Ohio. 

Simulation of short circuit faults on the synchronous alternator 
of a wind turbine generator, without resort to the classical assump- 
tions generally made for that analysis, indicates that maximum clear- 
ing times for the system tied to an infinite bus are longer than the 
typical clearing times for equivalent capacity conventional machines . 
Also, maximum clearing times are independent of tower shadow and wind 
shear. Variation of circuit conditions produce the modifications in 
the transient response predicted by analysis. 

The model and simulation developed for this study of transient 
stability are applicable for more extensive investigations of power 
systems involving wind turbine generators, and will continue to be 
useful until the extent that standard power system technology ap- 
plies directly to wind turbine power systems is fully determined. 
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ABSTRACT 


As a consequence of the world energy supply problems large kilo- 
watt and megawatt capacity wind turbine generators are being designed 
and built to take advantage of the energy in the wind. These appara- 
tus are equipped with synchronous alternators for direct synchroniza- 
tion with utility networks. Because a wind turbine generator inher- 
ently has unique power system features - such as a large moment of 
inertia relative to its power capacity and a continuously varying in- 
put torque - its electrical behavior must be investigated either by 
testing or by analysis. The transient response of the system to elec- 
trical faults is one mode of the system behavior that must be studied 
analytically rather than experimentally. 

In order to obtain a measure of its response to short circuits 
a large horizontal axis wind turbine generator is modeled and its per- 
formance is simulated on a digital computer.- The baseline model for 
the simulation is the DOE/NASA 100 kW wind turbine generator located 
near Sandusky, Ohio. 

Simulation of short circuit faults on the synchronous alternator 
of a wind turbine generator, without resort to the classical assump- 
tions generally made for that analysis, indicates that max im u m clear- 
ing times for the system tied to an infinite bus are longer than the 
typical clearing times for equivalent capacity conventional machines. 
Also, maximum clearing times are independent of tower shadow and wind 
shear. Variation of circuit conditions produce the modifications in 
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the transient response predicted by analysis. 

The model and simulation developed for this study of transient 
stability are applicable for more extensive investigations of power 
systems involving wind turbine generators, and will continue to be 
useful until the extent that standard power system technology ap- 
plies directly to wind turbine ■ power systems is fully determined. 
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Chapter I 


INTRODUCTION 

As the Department of Energy (DOE) of the United States expands 
its efforts to make wind turbine generators a viable and appreciable 
part of the national power network, it becomes necessary that power 
engineers gain a complete understanding of the electrical behavior 
of wind turbine generators-. The depletion of established sources 
of fossil fuels and the consequent price increases for those energy 
sources have turned attention to the energy available in sunlight 
and wind. 

1.1 Feasibility of Large Wind Turbine Generators 

There has been, geographically,- a widespread effort to produce 
commercially economical quantities of : electrical power with wind 
turbines.- Heretofore, all such efforts have been abandoned because 
of the relatively high cost of producing the electrical power from 
wind energy. 

The largest wind powered electric system that has been built 
to date was the Smith-Putnam machine, built in the early 1940’s in 
Vermont. This machine produced 1.25 MW of a.c. power (1) . The 
Danish government helped to develop and operate 200 kW experimental 
wind powered systems after World War II (2) . The British government 
built a 100 kW wind turbine generator in the Orkney Islands in 
1950 (3) . The French built and operated several large wind powered 
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powered electric generators in the period from 1958 to 1966 (4) . Some 
of the most advanced wind turbines that have yet been built were con- 
structed and operated in Germany during the period from 1957 to 1968 
(5). 

In every case the cost of the wind powered electricity was con- 
siderably in excess of the cost of conventional fossil fuel electric 
generation. As a consequence, the projects were abandoned. Despite 
their rejection what these early efforts to generate electricity from 
the wind energy proved was the technical feasibility of producing 
electrical power in this manner. 

1.2 U.S. Government Wind Energy Program 

The Federal Government has undertaken an accelerated program with 
the objective of stimulating the development of wind energy conversion 
systems capable of producing a significant amount of U.S. energy needs 
by the year 2000. To accomplish this, a series of experimental wind 
turbines from 100 kW to several megawatts rated capacity will be de- 
veloped over the next decade. It is anticipated that these develop- 
ments will be supported by an extensive research and technology effort, 
leading eventually to the installation of 10 to 100 MW multi-unit 
demonstration systems. 

The first large scale wind turbine generator was designed and 
built for the Department of Energy by the National Aeronautics and 
Space Administration (NASA) at the NASA/Lewis Research Facility at 
Plum Brook, Ohio, in 1976. One of the specific objectives of the 
NASA/Lewis wind energy program was to develop cost effective wind 
conversion systems that are compatible with user applications. The 
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wind power generator at Plum Brook is used as the baseline model for 
this research study. 

1.3 Power System Stability 

An- important -consideration of a practical power system is stabil- 
ity. ’ For an alternating voltage power system stability is that prop- 
erty which makes it po.ssible to- maintain synchronism between all in- 
terconnected machines for both steady state and transient conditions. 
The most important type of power system disturbance in transient sta- 
bility studies is a fault applied and subsequently cleared (6) . Sys- 
tem disturbances from faults may upset the balance between input and 
output - power of one- or more machines of a synchronous system. If the 
unbalance is. too great or too "prolonged, synchronism will be lost. 

The most effective method of improving stability is to reduce the 
duration of faults 'with the use of high speed relays and circuit 
breakers. This study of the effects of three phase short circuits 
- on a wind powered' synchronous alternator and- how some of the system 
parameters govern these effects has been made in order to predict the 
requirements- of fault correction relays and breakers for wind turbine 
power systems. 

1.4 Available Literature 

Although there exists an appreciable catalog of literature deal- 
ing with wind power (7, 8), there is little documentation of analyt- 
ical and experimental data related to the power system stability of 
wind power generators which have been built and operated in the past. 

t 

This sparse documented ^experience is the result of several conditions; 
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lack of funding to provide extensive analysis and data taking, private 
proprietary projects, lack of support for mathematical analysis, dis- 
posal of records when projects were discontinued, failure to catalog 
documentation. Further, the computational facilities and analytical 
techniques available today were not available to earlier projects. In 
any case, there is a scarcity of literature describing the electrical 
considerations and network experience of wind powered generators. 

With the advent of the current national program to develop wind 
turbines, analysis, of the electrical aspects of these systems has 
been undertaken. Hwang and Gilbert (9, 10) have studied synchroniza- 
tion of a wind turbine generator with an infinite bus. Johnson and 
Smith (11) have investigated the dynamics of wind powered generators 
on utility networks. Work has been done by the General Electric 
Company as part of an extensive analytical study for ERDA on wind' tur- 
bine system application.- Pantalone has studied effects of large in- 
terconnected wind generators on power systems (12) .- Hwang has also 

t 

studied the use of excitation control and reactance to help stabilize- 
wind turbine generators (13) . 

1.5 Distinctive Wind Turbine Generator Features 

There are distinctive features of wind turbine generators that 
make them different enough from conventional power systems to warrant 
an analytical study to determine the effect of wind turbine charac- 
teristics on power system stability. Among the characteristic fea- 
tures of wind turbine generators the two which have the most signif- 
icant effect on the power system performance are (1) high rotational 
inertia of the wind turbine, and (2) the variability of the turbine 
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excitation. The first characteristic is the product of the neces- 
sarily large area that the turbine blades must sweep in order to cap- 
ture appreciable energy from the wind. The second characteristic is 
the result of the inherent variability of the wind as well as modifi- 
cations to the wind stream that result from the large system struc- 
ture. 

This study of- the transient response of a wind turbine generator 
to - three phase faults considers especially these two characteristics 
and the concomitant features- that these characteristics produce. 
Variations of the baseline wind turbine model have been used to 
examine the effects of three phase faults on the performance of the 
electric power system. 

1.6 Summary 

The chapters which follow discuss the theory applied in this 
analytical study. Chapter II expresses the electrical theory in- 
volved in first swing transient stability analysis of synchronous 
machinery. Chapter III presents the model used for studying the be- 
havior of the wind turbine power system. Included in Chapter III is 
a discussion of the computer simulation. In Chapter IV the results 
of the analysis are presented and explained. The last section dis- 
cusses the conclusions derived from the analysis and the recommenda- 
tions for additional study of wind turbine power system stability. 
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ANALYSIS OF THE RESEARCH PROBLEM 

The transient stability analysis of a wind turbine generator can 
be approached with the standard analytical tools used to investigate 
alternating current electric power systems. The wind turbine genera- 
tor as any other alternating current power source, consists of a prime 
mover, a transmission, and an alternator with an excitation control for 
the alternator and a prime mover control to adjust the input torque. 
However, significant differences exist between the wind turbine gener- 
ator and its operating environment and the conventional power system 
and its environment. 

2,1 Research Problem 

The research problem is to determine the extent to which these 
differences affect the fault response of the wind turbine generator 
power system. 

Two distinctive features of the wind power system are considered 
in this analysis: 

(1) the turbine (prime mover) has much more rotary inertia than 

any other rotating component of the system, in fact, much 
more inertia than the rotating elements of conventional 
power systems of equivalently rated output power; 

(2) wind, the source of excitation (torque input) , is continu- 

ally varying. 
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A wind turbine generator is an electromechanical system, and any 
system disturbance produces a transient response. If any bounded in- 
put disturbance to the system produces a bounded output,- the system 
is stable. Power system stability is that property of a power sys- 
tem that enables it to remain in operating equilibrium (synchronism) 
during and following transient faults and circuit changes. 

The power system stability problem becomes very complex when a 
network of machines and loads is considered. A large set of non- 
linear coupled equations of the form 

x - f(x,u,t) (2.1) 

must be used to describe the system performance,. 

2.2 Problem of the Single Machine on an Infinite Bus 

' Before undertaking network problems, an understanding of the 
performance of a single machine must be developed. For this purpose 
one synchronous alternator tied electrically to an infinite feus is 
investigated. Figure 2.1 is a one- line diagram of a single alter- 
nator tied to an infinite bus through a transmission line. x £ Is 
the lumped reactance between the alternator voltage source and the 
infinite bus. Line resistance is assumed to be negligible. 

2.2.1 Swing equation. - Classically the motion of the alter- 
nator is described by the "swing equation’ 1 which relates the rotary 
inertia torque of the rotating shaft to the net mechanical and elec- 
trical torque on the shaft. 

I§. = T (2.2) 

For convenience the angle 8 is measured with respect to the 
synchronously rotating reference of the infinite bus. 




/ W5TT' 


Infinite Bus 



Figure 2.1- One machine tied to an infinite bus 
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6 = 6 + f a) Q t:' 

« tt 

e « s 

ti 

16 = T (2.3) 

The variation of 6 as a function, of time is indicative of the 
transient stability of the machine represented. In practice, the 
first swing of this function is usually the most significant in 
assessing the stability of the machine. If the rotor angle increases 
indefinitely, the machine will pull out of synchronism and stability 
will be lost. 

2.2.2 Classical assumptions. - The transient stability of one 
synchronous alternator tied to an infinite bus is a classical prob- 
lem. The solution is based customarily on the following assumptions 
( 6 ): 

(1) The mechanical input torque remains constant during the en- 
tire duration of the swing curve. 

(2) Damping or asynchronous power is negligible, 

(3) Synchronous power may be calculated from a steady state solu- 
tion of the network to which the machine is connected. 

(4) The machine can be represented by a constant reactance (di- 
rect axis transient reactance) in series with a constant voltage 
source (voltage behind transient reactance) . 

(5) The mechanical angles of the machine rotor coincides with 
the electrical phase of the voltage behind transient reactance. 

2.2.3 Applicability of assumptions to research problem . - For 
this research effort, which considers the performance of a wind tur- 
bine generator in response to an electrical fault, none of the 
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"classical assumptions" has been made. However, the applicability of 
the assumptions to the wind turbine generator problem is discussed. 

Assumption 1 is usually made because the input to the alternator 
is controlled by a speed governor. When a fault occurs at the output 
of the alternator, the output is abruptly changed, but the input re- 
mains unchanged. The alternator will vary about the synchronous 
speed, but until synchronism is lost, the speed change is very small. 
Usually the input governor will not respond until there has been a 
one percent speed change, and even then there is usually a time lag 
in the response. Accordingly, in the classical problem it is rea- 
sonable to assume that the mechanical input torque is constant dur- 
ing the duration of the swing curve. 

In the case of the wind turbine generator, however, the- torque 
input to the alternator is not controlled by a speed governor. The 
torque input to the alternator is a function of the power limit of 
the alternator. If the wind speed is below the system rated value, 
the rotor blades are set at a fixed pitch angle, and the torque input 
remains constant during a fault. If the wind speed is above the 
system rated value, the input torque is limited to its rated value 
by a power error signal controlling the rotor blade pitch. Then the 
input torque response to the output change resulting from the fault 
is limited only by the power sensor response time. 

Assumption 2 is usually an acceptable simplification in the 
classical model because damping or asynchronous power is small rela- 
tive to the synchronous power. The electrical output of a synchro- 
nous alternator consists of a synchronous part, dependent on the 
angular position of the machine, and an asynchronous part depending 
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on the angular speed of the machine. 

P e = P d <6) + P(6) (2.4) 

Although the asynchronous power may be small, its effect on tran- 
sient stability of wind power generators may not be negligible. As 
may be seen in plots of angle-time variation for a wind turbine gen- 
erator (Chapter. n£)., there is considerable variation of 6. The 
damping of the rotor swing is the result of the interaction between 
the airgap flux and the rotor windings, particularly the damper wind- 
ing, during slip ((d6/dt) ^ 0). 

Because the oscillation period of the alternator is relatively 
long in comparison to time constants of the network, steady-state 
conditions are usually assumed (Assumption 3) . In this research study 
two different network models are considered for comparison. 

The use of a fixed voltage behind the transient reactance for 
the alternator model (Assumption 4) is usually acceptable when the 
only electrical transient which must be considered is the transient 
component (rather than both the transient and subtransient components) 
and its time constant is longer than the period of mechanical oscil- 
lation. The transient component and its time constant are not neces- 
sarily longer than the mechanical oscillation period in the wind tur- 
bine generator with its large rotor inertia. The prolongation of 
the first swing due to the large rotor inertia and the availability 
of fast response voltage regulators mean that the assumption of 
fixed voltage during the transient is not warranted. Further, the- 
theoretical model of fixed voltage behind transient reactance is de- 
rived for the case of zero transient saliency or a round rotor ma- 
chine. To avoid any unwarranted simplification the alternator model 
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for this research is the complete set of Park's equations. The use 
of this model makes Assumption 5 unnecessary. 

The set of equations describing the model for this research was 
solved on a large digital computer. This facility made it unneces- 
sary and undesirable to make any assumptions merely to simplify the 
mathematical solution. 

Analyses of synchronous machines are simplified if the effects 
of saturation, hysteresis, and eddy currents in the iron are ne- 
glected. The importance of these neglected effects varies with the 
problem and the operating conditions. They are relatively unim- 
portant in the calculation of initial symmetrical rms short circuit 
currents of fundamental frequency, especially if the fault is sepa- 
rate from the synchronous machine by transformers and transmission 
lines (14) . 

2.2.4 Critical clearing time . - When a three phase dead short 
is placed across the output terminals of an alternator, the output 
power is reduced to zero immediately. There is then an accelerating 
torque applied to the alternator rotor because the input power has 
not been reduced. The alternator rotor continues to accelerate as 
long as this unbalance remains-. The displacement of the power angle 
with respect to the infinite bus increases and synchronism is lost. 

If the short is removed before synchronism is lost, the accel- 
erating power becomes negative with the resumption of the electrical 
load. If the decelerating power of the load is sufficient, the rate 
of change of the alternator power angle will become zero and reverse 
the angle change before synchronism is lost. For a given initial 
load, the maximum power angle at which the short circuit must be 
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cleared in- order to .prevent loss of synchronism is the critical clear- 
ing angle. Because power angle is ordinarily, iiotr known explicitly, 
critical clearing time, the time between application of the fault and 
clearing of the fault, is used in stability studies. The relays and 
breakers that are used to clear faults operate as timing devices 
rather than power . angle devices. 

2.3 Factors Affecting Clearing Time 

2. 3.1. Effect of inertia on clearing time . - Two synchronous 
generator design parameters which affect the critical clearing time 
for faults are the rotary inertia and the transient reactance. In- 
tuitively it is rational that an increase in inertia increases the 
time required to reach the critical angle. Analytically, from the 
classical swing equation it can be shown (6) that the time required 
for the power angle to change from 6^ to ^2 as a resu -*- t an ac “ 
celerating torque T is 



Usually no- literal solution for this equation exists. The 
equation indicates, however, the dependence of the clearing time on 
the square root of the inertia as well as on the accelerating ' 
torque. 

2.3^2 Effect of transient reactance on clearing time . - The 
effect of the transient reactance, x^, on the clearing time can be 
seen from the relationship that the accelerating torque in the swing 
equation is the difference between the input shaft torque and the 
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electrical torque. 


T = T sh - T e (2.6) 

For a synchronous machine with both a field and an amort is seur 
winding the fundamental frequency component of torque following a 
three-phase short is given by Concordia (15) as 


, -t/T. 

T (fund) = v 2 e 
e t 





+ — sin t 
*d 

(2.7) 


where . 


2x'*jx” 

Ta = + xj) (2 - 8) 

The first term within the brackets of equation (2.7) is the sub- 
transient component with a very short time constant relative to the 
transient component, the second term within the brackets. The first 
term can be neglected for the comparison being made. In the second 
term, x^ is ordinarily four to ten times the magnitude of x^. 
Therefore, the magnitude of the torque T decreases with increasing 
x^. The unbalanced accelerating torque thus increases with, increas- 
ing x^. As a consequence, the clearing time decreases with increas- 
ing reactance between the machine terminals and the short circuit. 
This conclusion follows from the fact that a symmetrical three phase 
reactance x e between the machine terminals and the short circuit 
..becomes part of the armature circuit and is added to x^, x^, x^j , x^, 
Xq> and x|j. 


2.3.3 Synchronous damping . - Two of- the functions of damper 
windings in a synchronous alternator are (1) to suppress hunting, 
and (2) to damp oscillations that result from faults. Damper action 
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in. the synchronous machine is explained with induction motor theory. 
In an induction machine the torque at small slips is almost directly 
proportional to the slip. Therefore, in calculation of synchronous 
damping the damping torque is assumed proportional to the slip. How- 
ever, because of the difference between the direct and quadrature 
axes, the constant of proportionality depends upon the angular posi- 
tion of the rotor, the power angle (16). An equation for the damping 
power has been developed by Park (17) and modified by Dahl (18) . 

The relationship is based on the following assumptions: 

(1) No resistance in either the armature or field circuit. 

(2) Small slip; where slip s -is .per unit slip-. 

(3) Damping caused by only one set of damper windings. 


o (x d " X P 
E 2 so) 1- TM 

(x e + xp 2 


sin 2 6 + 


(x 1 - x") 

q Q T »i 

(x„ + x*)2 W 


The derivative of this damping power with respect to <5 is 


f d 

d6 


-xp (x q - ^j) 

SW (x + x') 2 do ( X + x' )2 

jed 7 N e q' 

» 

(-4.2xl0~^ sin 26)6 for x e 

(-5.9x10“^ sin 26)6 for x g 


T" sin 26 

qo 

= 0.4 pu 
= 0.009 pu 


(2.10) 


Equation (2.9) shows the inverse relationship between damping power 
and external reactance which is one of the parameters examined in 
this research. Equation (2.10) has been evaluated for the synchro- 
nous alternator modeled in this study. The negative coefficients for 
practical values of the power angle 6 show that the damping power 
decreases with increasing load. 
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2.4 Wind Phenomena with- Respect to Power Output 

The wind, the normal driving energy source for a wind turbine 
system is variable in speed and direction. This variability is very 
much evident in the resultant drive train torque as well as the forces 
imposed on the structure of the wind turbine system.- Thus, even when 
the system is operating- with a fixed load', the input to the system 

* v 

is varying. Steady-state operation of the system is an infrequent 
mode and must be defined with reference to the input torque varia- 
tion. 

Several phenomena affect the velocity of the wind relative to 
the wind turbine, and, consequently, affect the torque developed in 
the drive train. These phenomena include wind shear, tower shadow, 
and blade flapping.. 

2.4.1 Tower shadow . - The significant wind variation in the 
technology of wind turbine generators is that caused by tower . 
shadow. 

Most of the large wind turbines that have been built have had 
the rotor downwind. The principal reasons for this configuration 
are that (1) the flapping blade is less likely to strike the tower 
on which the rotor is placed, (2) the rotor can be closer to the tower 
and less balancing mass is necessary to keep the center of gravity 
over the base of the tower,, and (3) if allowed to yaw freely, the 
rotor will seek a downwind location as a stable position. As a con- 
sequence of the rotor being downwind of the tower on which it is 
mounted, the tower is between the wind- and the rotor blades as they 
rotate. Even though the tower may be an open truss structure, it 
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provides considerable blockage to a free- stream of air. The impedance 
to the wind' and the resulting torque reduction in the turbine drive 
each time a blade passes behind the tower with- respect to the wind 
is termed tower shadow. 

The wind speed reduction and corresponding torque reduction that 
are the result of tower shadow can be as much as 35 percent with a 
truss tower (19) . The torque reduction occurs n times per rotor 
revolution for n-bladed rotor.* Thus, tower shadow produces an ever 

f 

present torque pulse n times per rotor revolution* 

2*4*2 Wind shear * - Wind shear is the result of ground surface 
friction on the wind* The wind speed normally increases with height 
above the ground. The speed variation is expressed by the relation- 
ship 

V » " V » r ef (h/ W k (2 ' 1U 

where the exponent k ( is a function of the reference surface rough- 
ness. 

As a result of the wind shear the torque produced by each blade 
varies as the turbine rotates. Because of the variation of the wind 
speed with height, there is a periodically varying shaft torque as 
the blades move through the wind field. For the wind turbine system 
studied a value of k = 0.22 is used. This is a k value repre- 
sentative of open country acceptable for installation of large wind 
turbines (20) . The torque variation resulting from the wind shear is 
not large (21) . The frequency of the variation is twice the turbine 
rotation frequency, however, and reinforces the tower shadow effect. 
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2.4.3 Blade flexibility . - The rotor blades of large wind tur- 
bines are relatively flexible structures. The response of these 
structures to aerodynamics and inertial loading results in blade 
motion, particularly flapping, which alters the torque developed. 
Study has indicated that these effects on rotor' shaft torque are 
not significant (21) . These effects are not considered in this 
study. 

2.5 Coupling of Tower Shadow Excitation and Output 

The tower shadow excitation has very important implications 'to 
the wind turbine generator designer and operator.. The entire sys- 
tem structure must be strong enough to withstand the loading which 
results from the pulsing excitation. The frequency of the excita- 
tion must be kept sufficiently different from any structural natural 
frequencies to avoid resonances. For the power engineers the con- 
cern is the extent to which tower shadow excitation is coupled to the 
system power output. 

2.5.1 Analysis of reduced model . - The existence of a twice per 
turbine rotor revolution torque pulse is inherent in a downwind rotor 
wind turbine with a two-bladed rotor. In order to reduce the cou- 
pling of that twice' per revolution excitation into the alternator 
power output the stiffness of the transmission drive can be reduced. 
The variation of the mode frequencies of the transmission drive train 
as a function of drive train stiffnesses can be examined with the 
simplified model shown in figure 2.2.. The figure represents the 
rotary inertias of the turbine rotor and the alternator coupled by a 
stiffness equivalent to all the drive train stiffnesses . The alter— 
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nator is represented by the second spring with a stiffness equal to 
the transient synchronous torque coefficient. The rotary motion of 
the transmission drive is given by equations (2.12) and (2.13). 

X 1®1 + K d< 9 l ~ e 2> = T in (2.12) 
I 2 0 2 ~ K d( 0 l ” e 2^ + k A 0 2 = i^sli (2.13) 


The eigenvalue frequencies of the model are given by equa- 
tion (2.14) and are plotted as a function of drive train stiffness 
in figure 2,3. 


9 _ K d + K A . 

Cij = 4* 

2I 0 


21 , 


(K d + K A ) [ K d (K d - K A ) ] 

.2 :zi,l 2 

4I 0 • 1 2 41, 


Sl/2 


(2.14) 


The magnitude of these frequencies are carefully considered by wind 
turbine structural designers in order to avoid a resonance from the 
tower shadow excitation of the input torque. In this study these 
resonant frequencies are of interest because the step action of a 

three-phase short circuit will excite these frequencies. 

Another, and a closer, approximation of these natural fre- 
quencies of the drive train and alternator can be obtained from an 
eigenvalue analysis of the closed-loop system using linearized ver- 
sions of the system equations. A linearized set of these equations 
is included in appendix B of this study. 


2.5.2 Effect of tower shadow on swing curve . - The tower shadow 
pulses that modulate the input torque of a downwind rotor wind tur- 
bine generator can affect the clearing time for faults across the 
output of the alternator. Although the excitation is unidirectional 
and approximates a half cosine wave in its variation, the alternator 
responds by swinging continuously, almost sinusoidally. (This re- 



Mode Frequency 1 (Hz) 



Mode Frequency 2 (Hz) 
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sponse is illustrated in Chapter lV'- where the results of the system 
simulation are described.) The rotor is continuously accelerating 
and decelerating. The clearing time for a three-phase short circuit 
depends upon the rotor acceleration due to the power pulsations plus 
the acceleration due to the unbalanced torque from the fault. 

2.5.3 Effect of alternator excitation system on swing curve . - 
The effect of the alternator excitation system on the severity of 
the first swing is relatively small, except for faults with long 
clearing times (22). Avery fast, high-response, excitation system 
will usually reduce the first swing by only a few degrees, or will 
increase the generator transient stability power limit by a few per- 
cent. The large inertia constant H of a wind turbine alternator 
makes for longer than usual clearing time. A representation of the 
alternator excitation system is warranted for this study. 

2.6 Summary 

The type of power system disturbance which is most important 
in stability studies is a fault applied and subsequently cleared 
(6). Faults are usually not permitted to remain on power systems 
any longer than is required for relays to detect their presence and 
open the appropriate circuit breakers. 

The classical approach to the transient stability analysis of 
a faulted synchronous alternator tied to an infinite bus involves 
several simplifying assumptions. With the availability of modern 
digital computers for the solution of sets of nonlinear differential 
equations one of the reasons for making the assumptions is elimi- 
nated. At this early stage of investigations into the developing 
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technology of wind turbine generators,- it is preferable to represent 
the system characteristics more closely. Concordia points out (23) 
that when the effect of delayed fault clearing is to be studied, a 
detailed alternator representation is desirable even though only the 
first swing transients of the alternator rotor are considered. If 
the simple concept of generator voltage behind transient reactance 
is used, the effect of clearing time is smaller than in reality. It 
follows, that critical clearing times determined by a simplified sta- 
bility survey are likely to be longer than in reality. 



Chapter III 


SIMULATION MODEL 

The model used in this study is that of a typical large hori- 
zontal axis wind turbine generator system. The essential elements 
of the wind turbine system model are shown in the block diagram in 
figure 3.1. This is the configuration of the DOE/NASA 100 and 200 KW 
wind turbine generators already constructed as well as the larger 
2 MW systems in development for the DOE. 

A representative model consists of a large inertia turbine which 
converts wind energy to mechanical energy by applying torque to a low 
speed (20-40 r/min) shaft. A speed increaser converts the low speed, 
high torque energy to a high speed (1800-3600 r/min) , low torque en- 
ergy required to drive a conventional synchronous alternator. This 
transmission is illustrated in figure 3.2. 

For this study the model has a rigid turbine developing torque 
that is transmitted by a drive shaft through a 1:45 speed changing 
gear box. The reduced torque is transmitted to a synchronous gen- 
erator. The model is detailed in the following sections. The 
def initions-df the terms in all equations are given in appendix A. 

3.1 Turbine 

The turbine model is a two-bladed propeller. The magnitudes of 
the useful torque developed by the turbine 1 .ar'e' calculated with a 
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Figure 3.2. - Energy transmission through wind turbine generator. 
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digital computer program developed by Wilson and Lissaman (24) and 
modified for this research program. The turbine torque computer pro- 
gram uses a modified blade element theory and includes corrections for 
tip losses. The usable wind energy depends on wind velocity, turbine 
rotary speed, blade pitch angle, airfoil characteristics, and tur- 
bine geometry.- The turbine geometry consists of the characteristics: 
blade length, number of blades, coning angle, and radial distributions 
of chord and twist angle. The turbine model computer program calcu- 
lates shaft power, shaft torque, and thrust on the blades. 

For a given blade geometry the torque iisi generated as a func- 
tion of the wind speed perpendicular to the plane of the turbine, the 
turbine rotational speed, and the blade pitch angle. The relation 
for torque is 

T in - f <V V B l> «•« 

Torque-pitch angle characteristics calculated by the program for sev- 
eral wind speeds at 40 r/min are shown in figure 3.3. 

3.2 Drive Train 

The drive train elements, that is, the low speed shaft', the gear 
box, the high speed shaft, and the couplings,- have been modeled as the 
pi-network indicated in figure 3.4.- The rotational inertias of all 
elements on the low speed side of the gear box, including the gear 
box, have been combined into one inertia, 1^, separated from the com- 
bined inertia, I^, of all components on the high speed side of the gear 
box, including the generator rotor, by the composite stiffness, K^, of 
all elements between the turbine and the alternator. Transmission 
train lossed are accounted for by the damping factors and B^. 






0i &z 



Figure 3. 4- Equivalent Pi-network of Drive Train 
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The simplification of the drive train to a pi-network is justi- 
fied on the basis that 

(1) the rotational inertia of the turbine is much greater 

(15 times) than the combined rotational inertias of all 
other elements of the drive train, and therefore, it alone 
essentially determines I. , 

(2) the alternator rotor provides the largest rotational inertia 

on the high speed side of the gear box, and essentially 
determines ly , and 

(3) the large low speed drive shaft and the gear box are the 

softest springs in the drive train and essentially deter- 
mine its effective stiffness. 

The stiffness of the drive train determines the amount of cou- 
pling between the turbine and the alternator. The extent of this 
coupling becomes important to the network performance of the wind 
turbine generator because -of the tower shadow effect on the - input 
torque excitation (21) . Composing all .the drive train stiffnesses 
into one spring facilitates treating the coupling parametrically while 
permitting retention of the lowest oscillatory mode of the drive train 

' r 

which is in a range most likely to excite power network modes. Since 
the stiffnesses and inertias of elements other than the turbine, the 
low speed shaft, and the alternator rotor are not readily available 
from measurements or calculation, they are estimated. Little ac- 
curacy is lost, therefore, in resorting to an equivalent pi-network 
for the drive train. The numerical values of the baseline model 
drive train parameters are listed in Table 3.1. 



31 


Table 3.1 


STIFFNESS AND MASS MOMENT OF INERTIA OF 
BASELINE DRIVE TRAIN COMPONENTS 



Mass moment 

Torsional 


of inertia, 

stiffness 


5 > 

N-m-sec 

N-m/rad 

Turbine 

127,735 

"Rigid" 

Low speed shaft 

835 

2.26xl0 7 

Gear box 

281 

1.88xl0 7 

High speed shaft a 

b 992 

11.54xl0 7 

Alternator rotor a 

3,887 

96.40x10' 

Pulleys and belts 3 

3,356 

7.28xl0 7 


a Ref erred to turbine speed. 
^Includes disk brake. 


The state space model of the equations of motion of the drive 
train is 




(3.2) 


(3.3) 



(3.4) 


3.3 Alternator 

The electrical generator modeled is a synchronous alternator de- 
scribed the classical Park’s equations (17) » modified by Olive (25) 
to eliminate any explicit expression of inductance. • The alternator 
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equations reflect the assumption that the synchronous generator fre- 
quency is the same as the bus frequency as long as synchronism is 
maintained (26) . No assumptions are made with respect to the relative 
magnitudes of the "transformer voltages" (dX^/dt, dX^/dt) and the 
"speed voltages" (wX d , (oXq) . The following equations describe the 
alternator. 



e" - X . 
i = —9 1 

d X d 

(3.5) 


e" + X 

± a 

JL It 

(3.6) 


q -x d 


dXj 

U _ 

dt 

[v d + Vd ‘ (e d + x ,V“ ] “o 

(3.7) 



[v + r i + (x.i, - e i - e 0 )w]w n 
1 q a q d d ql q2 0 
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_ q q q q d 

e d ■ x" 
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< x d ■ X P 


< x d “ x d } 


"ql (x' - xj) q (x* - xj) q 


e_o - 


'q2 




X P*d 


+ x"e* 


x'e" 


-x 


< x d ~ x d> 

< x d - X P 


d q d q 


de* 

= 

dt 



de a (x d “ X P 
dt " " ' X P 


(3.8) 


(3.9) 


( 3 . 10 ) 


(3.11) 


(3.12) 


(3.13) 
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The electromagnetic torque is 


T = e"i, + e"i - i i (x" - x") 
e ad qq a q a q 


(3.14) 


m = TI / poles U 737.8 kvabase 

L V 2 A » 6 j e 


lb -ft 


737.8 kvabase' 


(0.1383)T e n-m 


The machine reactances appearing in the alternator equations 
include the output transformer and the tie line reactances. There- 
fore, the terms v. and v in equations (3.7) and (3.8) are the 
— q 

infinite bus direct and quadrature axes voltage components as shown 
in. figure 3.5 where v fc is the alternator terminal voltage. 


3.4 Alternator Excitation System 


The modeled excitation system tends to maintain a constant power 
factor load on the alternator as well as a constant output voltage 
(13). -Shown in block diagram form in figure 3.6 the excitation sys- 
tem is described by the following equations: 


. 4v t - V r.f - + Srfc - PF • A a) 


Av f = — Av t I exp - — 


+ 


where, as shown in figure 3.5, 



original pa gf, t<? 

op pooe a (S 



Figure 3.5-Terminal Voltages 
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Figure 3,6. - Constant power factor excitation system 
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td V d ^q X e 

(3.20) 

• - v + i.x 

(3.21) 

tq q d e 

■ e- * <r 

(3.22) 

v d « sin 6 

(3.23) 

v - cos 8 

q 

(3.24) 


The term (v^ ef - v fc ) in equation (3.17) represents the voltage devia- 
tion; the term ((p/v ) - PF » i ) represents the power factor devia- 

Z 2L 

tion. 

3.5 Turbine Blade Pitch Control 

Blade pitch control is used to regulate power output of the al- 
ternator. At wind speeds below the system rated wind speed, alterna- 
tor output is correspondingly below the rated alternator power. At 
these less than rated conditions the blades are kept at a fixed pitch - 
'•maximum power position." When the wind speed is above the rated value 
(Table 3.1)', the alternator is actively limited. Limitation is ef- 
fected by pitching the turbine blades. Output power is measured and 
compared with the reference value to generate an error signal that 
actuates the controller and the blade pitch servomechanism (fig. 3.7) . 
The equations describing the control action at and above rated wind 
conditions are listed in state space form. 
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Figure 3.7-Cbntroller and Blade Pitch Servo Models 
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3.6 Computer Simulation 


(3.25) 


The wind turbine generator model described in the preceding sec- 
tions provides a set of nonlinear, interrelated, differential equa- 
tions. This set of equations was solved numerically with a computer 
code developed for this study using CSMP - Continuous System Modeling 
Program (27) . Integration is performed using a fourth order Runge- 
Kutta integration method with a fixed step size of 0.002. Options are 
provided in the code for the inclusion of tower shadow effect, wind 
shear effect, and gusting of the wind. 

The CSMP INITIAL segment of the code includes a set of equations 
that establishes the initial conditions of the system for any selected 
power and power factor loading of the alternator. In the CSMP 
DYNAMIC segment of the program a three phase short circuit can be ap- 
plied at any time for any duration. Subsequent to the application of 
the fault the bus and alternator rotor speeds, which are identical 
prior to the occurrence of the fault, are computed independently. 

In cases for which a constant wind speed is assumed, the output 
transients result from a three phase fault short circuit applied- at 
time t = 0. In cases for which tower shadow effect is included, the 
transient resulting from the introduction of tower shadow is permitted 
to stabilize to a 'feteady state" prior to the application of the fault. 

A listing of the computer program used for this study is included 
as Appendix F of the dissertation. 
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3.7 Summary 

A block diagram of the complete wind turbine generator system 
model for this study is shown in figure 3.8. The model has been con- 
figured for simulation on a digital computer to provide a measure of 
the response of the alternator to three phase zero impedance short 
circuits . 

The numerical values of the system constants used in the simula- 
tion are listed in Appendix E. 




Figure 3.8“Simulation model block diagram 










Chapter IV 


RESULTS OF SIMULATION 

The response of the synchronized wind turbine alternator to a 
three phase, zero impedance, short circuit was simulated to examine . 
the nature of the swing curve and to determine maximum clearing time 

4 

required to maintain synchronism. A tabulation of the numerical re- 

* 

suits of the simulation are presented in Table 4.1 for the selected 
parameters and operating conditions examined in the study. 

4.1 Selected Operating Conditions and Parameters 

4.1.1 External reactance . - Two values of external reactances 
were selected for comparison. The lower value, x e = 0.009 pu, is - 
the tie line reactance of the Mod-0 wind turbine which was used as 
the baseline model for the research study. The second tie line re- 
actance, x e = 0.4 pu was selected as representative of a high reac- 
tance so that the two reactances spanned the probable range of values 
to be expected realistically. In both cases the line resistance was 
assumed negligible. 

4.1.2 System parameters . - Two parameters of the wind turbine 
system, the equivalent drive train stiffness and the total rotary in- 
ertia, were varied to observe the effects on first swing transients 
resulting from short circuit faults. Both of these parameters are 
structural features of concern to system designers as well as power 
engineers 
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TABLE 4.1 

• MAXIMUM CLEARING TIME (WITHOUT LOSS OF SYNCHRONISM) 


Tie-line 

Drive train 

Relative 

Number 

of 60 Hz cycles in 

reactance. 

relative 

rotary. 

maximum short 

at given 

pu 

stiffness 

inertia 1 ’ 


wind speed 

• 


• 

14 mph 

18 mph 

30 mph 






(wind speed*) 

0.009 

0.1 

1.0 

25 

15 

15 

* 

1.0 

1.0 

67 

37 

37 


10.0 . 

1.0 

67 

38 

38 


.1 

.5 

18 

12 

12 


1.0 

.5 

‘ 5.6 

28 

29 


10.0 

.5 

53 

29 

30 

0.4 

0.1 

1.0 

15 

7 

7 


1.0 

1.0 

45 

21 

21 


10.0 

1.0 

44 

23 

24 


.1 

.5 

11 

5 

5 


1.0 

.5 

33 

15 

16 


10.0 

.5 

32 

17 

18 


*Wind speed 

14 mph = 6.2 m/s 
18 mph = 8.0 m/s 
30 mph = 13.4 m/s 


Corresponding load 

0.4 pu, 0.8 PF 
0.8 pu, 0.8 PF 
0.8 pu, 0.8 PF 


4.1.3 Electrical loads . - Three loading conditions were se- 
lected for examination: a 0.4 pu power load at a wind' speed of 

14 mph, a 0.8 pu power load at a wind speed of 18 mph, and a 0.8 pu 
power load at a wind speed of 30 mph. All the loads have a 0.8 
lagging power factor. The 0.8 pu load at 18 mph wind speed repre- 
sents the rated loading at rated wind speed for the baseline model 
wind turbine. The 0.4 pu load is a 50 percent rated load and occurs 
in the fixed blade pitch control mode of the system operation. The 
0.8 pu load at 30 mph wind is the rated load in the variable pitch 
or power limiting control mode of the system. The rated load at rated 
wind speed is nominally at the interface of the two control modes . 
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4.2 Determination of Maximum Clearing Tims 

To determine the maximum clearing time for a particular set of 
operating conditions and system parameters a three phase short circuit, 
zero impedance, fault was applied to the system at the infinite bus 
terminals. The duration of the simulated fault, measured in terms 
of cycles of the 60 Hz frequency, was increased successively in in- 
tegral values until a power angle swing that exceeded 180° was ob- 
served. The short circuit transient preceding the occurrence of the 
slipped pole determined. .the maximum clearing time. 

4.3 Tower Shadow and Wind Shear 

The ever present tower shadow phenomenon in the downwind rotor 
wind turbine generator and the wind shear produce a twice per tur- 
bine rotor revolution torque excitation that is coupled to the alter- 
nator through the drive train transmission. (A two-bladed turbine 
rotor is implied.) The effect of this excitation on the baseline 
system power output. ' is shown in figures 4.1 and 4.2. 

A three phase short circuit fault applied to the system at a 
time when the alternator rotor acceleration is at a maximum produces 
the transient displayed in figures 4.3 and 4.4. The figure displays 
the maximum clearing time for the system operating at the given power 
conditions and experiencing tower shadow and wind shear. A short 
circuit duration 1 cycle longer causes loss of synchronism, and 
that result is also indicated in figures 4.3 and 4.4. 

A short circuit applied to the baseline model simulated without 
the tower shadow and wind shear effects produces the transient and 
maximum clearing time illustrated in figures 4.5 and 4.6. The com- 
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Figure 4.1. - Baseline model power angle oscillations from tower shadow and wind shear 







Figure 4,3. - Swing curve of baseline model with tower shadow and wind shear. 
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Rotor speed curves of basel-ine model with tower shadow and wind shear./ 
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Figure 4.6. - Rotor speed curves of baseline model without tower shadow and wind shear 
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parison is more easily seen in figures 4.7 to 4.10 in which the tran- 
sient curves of figures 4.3 to 4.6 are superimposed. The similarity 
of the two transients implies that for the study of the' effects of 
this fault on the wind turbine generator tower and wind shear can be 
neglected. This inference is drawn because the magnitude of both 
tower shadow and wind shear simulated is probably as high as each 
would ever be in a real wind turbine generator system. 

For the nominal turbine speed used for the model the twice per 
turbine rotation excitation produces a 1.33 Hz excitation frequency. 
This frequency is above any of the low frequency mode values indicated 
in Table 4.2 and does not excite other oscillation modes. 

TABLE 4.2 

COMPARISON OF SIMULATION AND EIGENANALYSIS FREQUENCIES 
Alternator stiffness, 2537 lb-ft/elec rad 1390 lb-ft/elec rad 


k a 

(x e 

= 0.009) 

(x e = 0.4) 


Drive train relative 






stiffness, K^ a 

0.1 

1.0 

10.0 

0.1 1.0 

10.0 




Low-frequency mode 


Simulation frequency, 
Hz b 

Eigenanalysis fre- 

0.3 

0.7 

0.7 

0.3 0.6 

0.8 

quency, Hz 

0.3 

0.7 

0.8 

0.3 0.6 

0.8 




High-frequency mode 


Simulation frequency, 
Hz b 

Eigenanalysis fre- 

5.1 

6.4 

14.0 

3.9 4.8 

13.3 

quency, Hz 

4.8 

6.0 

13.5 

3.6 5.2 

13.1 


a 

Stiffness relative to baseline model stiffness. 
bFrequency determined from digital computer output listing. 

The option to neglect the tower shadow and wind shear effects 
was selected for this research study. This simplification materially 
reduced computer time. If there is error in making this assumption. 
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With tower shadow, 40 cycles 
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Wind speed 30 mph 
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I 1.0 

X c 0.009 PU 
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Figure 4 . 10 , - Comparison of rotor speed curves of 'baseline model with and without 

TONER SHADOW AND WIND SHEAR SHOWING LOSS OF SYNCHRONISM. 
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it is in the direction to reduce the maximum clearing time; that is, 
the determination is conservative. 

4.4 Typical Swing Curve 

The typical result of the scheme of clearing time determination 
adopted for this study is illustrated in figure 4.5. In this ex- 
ample < a 37 cycle duration short circuit produces one big swing 
that rapidly is damped so that the power angle returns to its steady 
state value. A short circuit one cycle longer produces sufficient 
acceleration to the rotor to cause loss of synchronism. The same 
alternator rotor behavior is illustrated in figure 4.6 which shows 
the relative speed of the rotor. Upon application of the fault at 
time = 0, the rotor accelerates. At time 0.617 second (37 cycles) 
the fault is removed and the frequency difference between the al- 
ternator voltage and the bus voltage produces the ensuing speed 
oscillation that rapidly is damped. Increasing the short circuit 
duration to time = 0.633^ second (38 cycles) produces sufficient rotor 
acceleration to drive the alternator out of synchronism. 

The configuration of the swing curve presented in figure 4.5 and 
the damped oscillation of the rotor speed curve of figure 4.6 are 
representative of the transient responses produced by the application 
of a three-phase short to the wind turbine generator system output 
for different system parameter values and/or operating conditions. 

In each case for which synchronism is not lost, application of the 
fault produces a large positive swing of the power angle curve that 
lasts for 1 to 3 seconds before the angle returns to a damped oscil- 
lation about the final value. The generalization is indicative of the 
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order of magnitude of the duration of the transient. 

4.5 Oscillation Modes 

Two frequencies are evident in the swing and the ensuing damped 
portion of the transient. The same oscillation frequencies are 
recognizable- in the rotor speed curves. The higher frequency 
oscillation is rapidly damped; the lower frequency experiences less 
damping as the transient gradually subsides to the steady state 
value . 

Changing system parameter values and/or operating conditions 
changes the frequency and amplitude of the oscillation, but does not 
alter the transient function radically. 

The eigenvalue analysis of the simplified model depicted in fig- 
ure 3.4 provides good predictions of these two modes of oscillation 
* which occur in the more detailed, closed loop simulation output. 

The correspondence of the system simulation output frequencies de- 
termined from the digital computer printout and the eigenvalue anal- 
ysis frequencies is shown in table 4.2. 

The lower frequency mode is the effective natural frequency of 
the alternator. It is determined by the "spring" of the alternator 
synchronous torque and the rotary inertia of the rotating shaft. 

The higher frequency mode is produced by the spring-mass combination 
of the alternator rotor and the synchronous torque paralleled to the 
drive train spring.- 

4.6 Effect of Rotational Inertia 


The influence of the rotary inertia on the maximum clearing time 
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is shown in Table 4.3. Tabulated are the ratios of maximum clearing 
times for system conditions matching in all respects except rotary in- 
ertia.- The mean value for the 18 ratios examined is 1.3. The pre- 

1/2 

dieted ratio from the analysis (eq. (2.5)) is 2 . The result of 

the simulated effect of inertia on clearing times matches the ex- 
pectations from earlier analysis. 

4.-7 Effect of Power Load 

Clearing time is effectively inversely proportional to load . 

The dependence of- the duration of maximum clearing time on alterna- 
tor power, level is illustrated by the ratios in Table 4.4. The ra- 
tio of clearing times for a 0.4 pu load and an 0.8 pu load for 12 
otherwise matching systems is- calculated. The mean value of the 
12 ratios tabulated is 1.9, corresponding to the inverse of the load 
ratio. This result implies the generalization that, in spite of the 
nonlinear ities present in the wind torque characteristics, in the 
alternator model, in the control modes, and in the parameter selec- 
tions used, clearing time is effectively inversely proportional to 
load. The, implication is further substantiated by the result that 

I 

clearing times for an 0.8 pu load at 18 mph and at 30 mph wind speeds 
are substantially the same. 

This close correspondence in clearing time ratios exists even 
though the values of the partial derivatives which describe the wind 
turbine torque characteristics differ considerably at the two wind 
speeds investigated. The differential torque developed by the wind 
turbine can be expressed by equation (4.1). (See also Appendix D.) 



TABLE 4.3 


EFFECT OF ROTARY INERTIA ON MAXIMUM 


Tie-line 

Drive train 

Relative 

Number 

-of 60 Hz 

cycles 

reactance. 

relative 

rotary 

in maximum short at 

pu 

stiffness 

inertia 

given wind speeds 




14 mph 

18 mph 

30 mph 

0.009 

0.1 

1.0 

25 

15 

15 


.1 

.5 

18 

12 

12 


1.0 

1.0 

67 

37 

37 


1.0 

.5 

56 

28 

29 


10.0 

1.0 

67 

38 

38 


10.0 

.5 

53 

29 

30 

0.4 

0.1 

1.0 

15 

7 

7 


.1 

.5 

11 

5 

5 


1.0 

1.0 

45 

21 

21 


1.0 

.5 

33 

15 

16 


10.0 

1.0 

44 

23 

24 


10.0 

.5 

32 

17 

18 


*Wind speed Corresponding load 

14 mph =6.2 m/s 0.4 pu, 0.8 PF 
18 mph = 8.0 m/s 0.8 pu, 0.8 PF 
30 mph = 13.4 m/s 0.8 pu, 0.8 PF 


TIME 


Clearing times ratio for 2:1 
inertia ratio 

14 mph 18 mph 30 mph 

(wind speed*) 

1.4 1.2 1.2 

1.2 1.3 1.3 

1.3 1.3 1.3 

1.4 1.4 1.4 

1.4 1.4 1.3 

1.4 1.4 1.3 


Mean ratio = 1.3 
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0.009 


0.4 


TABLE 4.4 

EFFECT OF POWER LOAD, ON MAXIMUM CLEARING TIME 

Drive train Relative Number of 60 Hz Clearing times 

relative rotary cycles in maxi- ratio for 0.4/0. 8 

stiffness inertia mum short at load ratio 

given wind 
speed* 

14 mph 18 mph 


0.1 

1.0 

25 

15 

1.7 

1.0 

1.0 

67 

37 

1.8 

10.0 

1.0 

■67 

38 

1.8 

.1 

.5 

18 

12 

1.5 

1.0 

.5 

56 

28 

2.0 

10.0 

.5 

53 

29 

1.8 

0.1 

1.0 

15 

7 

2.1 

1.0 

1.0 

45 • 

21 

2.1 

10.0 

1.0 

44 

23 

1.9 

.1 

.5 

11 

' 5 

2.2 

1.0 

.5 

33 

15 

2.2 

10.1 

.5 

32 

17 

1.9 


*Wind speed 

14 mph =6.2 m/s 
18 mph =8.0 m/s 


Corresponding load 

0.4 pu, 0.8 PF 
0.8 pu, 0.8 PF 


Mean ratio =1.9 


AT 


in 


3T. 9T . 3T. 

T“ AV + A6 + AB, 

3V w .so 1 38 1 

w d,:, 1 


- K V AV w + A0 1 + K B Ae i 


(4.1) 


The variation of the values of the gain factors (partial derivatives) 
for the wind speeds 14, 18, and 30 mph are shown in Table 4.5.- A 
plot of the wind turbine torque characteristics are shown in fig- 
ure 4.11. 

The relationship between alternator load and clearing time is 
seen in the similarity of the transient swing curves and speed curves 



TURBINE TORQUE s> 
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TABLE 4.5 

GAIN VALUES OF DIFFERENTIAL INPUT TORQUE 
AT MODELED WIND SPEEDS 


Wind speed, 

K v , 

Kfl. 

Kg, 

mph 

(f t-lb/ft/sec) 

(f t-lb/rad/sec) 

(f t-lb/rad) 

14 

1771 

-4 152 

71 418 

18 

2441 

-4 834 

121 127 

30 

2946 

-20 382 

393 707 


for the three sample wind speeds applied to the same system il- 
lustrated in figures 4.12 and 4.1-3. The clearing time for the 
14 mph, 0.4 pu load is almost, twice as long as the 0.8 pu loads. 

The same comparison is made for the three sample loading conditions 
in figures 4.14 and 4.-15 for a tie-line reactance of 0.4 pu. For 
each value of x e the reduced damping associated with the lower 
frequency for the 0.4 pu power trace is evident.- The lower fre- 
quency is associated with the alternator. The more lightly loaded 
alternator is more lightly damped. This result is consistent with 
the theory expressed by equation (2.10). 

4.8 Effect of Drive Train Stiffness 

Reduced drive train stiffness has a "decoupling" effect. 

The effect of different drive train stiffness is illustrated in fig- 
ures 4.16 and 4.17, one set of curves for x e = 0.009 pu, and the 
other set for x e = 0.4 pu. In each set of curves the transient 
swings for the system with the baseline value of the drive train- stiff- 
ness and the system with the increased value of ten times the baseline 
stiffness are not very different. However, the systems with the re- 
duced drive train stiffness show the "decoupling" effect of the 
softer transmission line. In the case of each external reactance 
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Figure 4.17, - Swing curves for "soft" tie-line with selected drive train stiffnesses^ 
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the same severe three phase, zero impedance, short circuit that 
causes the large swings for the stiff er systems, produces in the 
soft systems a transient whose frequency components are damped in 
one- third the time as the frequency components of the stiffer sys- 
tems. 

The ''decoupling” that results from the soft transmission also 
has the effect of reducing the maximum clearing time over 60 percent. 
The large turbine rotary inertia is not as closely coupled as in the 
stiffer systems, and as a consequence of the reduced inertia the 
rotor accelerates faster. Synchronism is lost more quickly. This 
effect can be noted from a comparison of figures 4.18 and 4.19. 

4.9 Effect of Tie-Line Reactance 


The appreciable increase of the reactance between the fault and 
the alternator terminals from x e - 0.009 pu' to x g = 0.4 pu reduces 
the maximum clearing time by almost 50 percent for otherwise corre- 
sponding circuit conditions. For the 18 pairs of corresponding con- 
ditions tabulated in Table 4.6 the average of the ratios of clearing 
time for x e - 0.4 to clearing time for x e = 0.009 is 0.56. This 
reduction is in accord with the reduced synchronous torque of the al- 
ternator. The equation for the synchronous torque is given by 


e q e 2 


sync x! + x o 


cos 6 + 


e 2 (x r - 
e 2 ' x d 




(x' + x e ) (X q + x e ) 


cos 25 


(4.-2) 


The effect of increasing the reactance x g is to reduce the torque 
coefficient T • The torque coefficient is effectively the stiff- 
ness of the spring indicated in figure 2.2. Reducing the stiff- 

ness of the system enables the rotor to accelerate more rapidly, and 
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Figure 4,18, - Swing curves for "stiff" tie-line with selected drive train stiffnesses 

SHOWING LOSS OF SYNCHRONISM. 
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TABLE 4.6 


EFFECT OF TIE-LINE REACTANCE ON MAXIMUM CLEARING TIME 


X 

e 

Drive train 
relative 
stiffness 

Relative 

rotary 

inertia 

Number of 60 Hz cycles 
in maximum short at 
given wind speeds 

14 mph 18 mph 30 mph 

Clearing time ratio for x g 
increase 

14 mph 18 mph 30 mph 

(wind speed*) 

0.4 

.009 

0.1 

.1 

1.0 

1.0 

15 

25 

7 

15 

7 

15 

0.60 

0.47 

0.47 

0.4 

.009 

1.0 

1.0 

1.0 

1.0 

45 

67 

21 

37 

21 

37 

0.67 

0.57 

0.57 

0.4 

.009 

10.0 

10.0 

1.0 

1.0 

44 

67 

23 

38 

24 

38 

0.66 

0.60 

0.63 

0.4 

.009 

0.1 

.1 

0.5 

.5 

11 

18 

5 

12 

5 

12 

0.61 

0.42 

0.42 

0.4 

.009 

1.0 

1.0 

0.5 

.5 

33 

56 

15 

28 

16 

29 

0.59 

0.54 

0.55 

0.4 

.009 

10.0 

10.0 

0.5 

.5 

32 

53 

17 

29 

18 

30 

0.60 

0.59 

0.60 







Mean ratio 

- 0.56 

*Wind speed 

Corresponding load 





14 mph 
18 mph 
30 mph 

- 6.2 m/s 
= 8.0 m/s 
= 13.4 m/s 

0.4 pu, 
0.8 pu, 
0.8 pu. 

0.8 PF 
0.8 PF 
0.8 PF 
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consequently the maximum clearing time for the fault is reduced. 

The reduction of the alternator "stiffness" is indicated also 
in Table 4.2 by the reduced frequencies for the x e = 0.4 cases. 

This reduction is in accordance with equation (2.12) which relates 
the mode frequencies to the drive train and alternator stiffness fac- 
tors. 

The additional line reactance has reduced the transient stabil- 
ity of the power system so that for a given steady state load the 
system becomes unstable for a shorter fault duration. This behavior 
of the system is displayed in figures 4.20 and 4.21. 

Increasing the tie-line reactance x £ also reduces the damping 
power that is produced by the fault excitation. This reduction is 
consistent with the analytical expression for damping power 
.(eq. (2.9)). 

The changes in maximum clearing time and in the damping of the 
transient following a short circuit to two systems with different 
external reactances are depicted in figures 4.22 to 4.27, For a given 
steady state load the system becomes unstable for a shorter fault 
duration as figures 4.20 and 4.21 show. 

4.10 Effect of Neglecting "Transformer Voltages" 

As implied in Section 2.2.2 the classical assumptions for cal- 
culating the transient response of one synchronous alternator when 
faulted with a three-phase short neglect the so-called transformer 
voltages (dX^/dt, dX^/dt) on the basis of the assumption that they 
are small in comparison to the so-called speed voltages (ojX^, coX^) . 
This assumption was not made during the previous simulations. 
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Figure 4.20. - Swing curves for "stiff" tie-line at selected wind speeds showing loss 
of synchronism. 
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Figure 4.22. - Comparison of swing curves of baseline model with selected tie-line 

REACTANCES . 
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ACTUAL SPEED/RATED SPEED 


Wind speed 30 mph 

k d i.o 

I 1.0 


1 i 1 i » A 
* ill 11 A A 

II J ^ 1) ^Xp 0. PO, 21 CYCLES 

1 '| l| ‘ 

1 I i< 


~ X_ 0.009 pu, 37 cycles 


1.3* l.S 2.3 2*5 3.3 3.S ‘Ktf S-.S S.tf 

TIME SEC 

Figure 4.23. - Comparison of rotor speed curves of baseline model with selected tie-line 

REACTANCES. 
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Figure A, 25, - Comparison of rotor speed curves of "soft" drive train model with 

SELECTED TIE-LINE REACTANCES, 
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Figure 1 , 27 , - Comparison of rotor speed curves of "stiff" drive train model with 

SELECTED TIE-LINE REACTANCES, 
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In order to compare the transient response with and without the 
simplifying assumption, the research computer program was modified 
to assume the transformer voltages are negligible. When this assump- 
tion is made, the alternator equations (3.7) and (3.8) are modified to 
equation (4.3) and (4.4). - ' 


i d = 


’ r (v , - e") + x" (v - e") 

_ . A -i £. 3_3 31 

x"x" + r 2 
d q a 


(4,3) 


- 


X d (V d ' e P ~ r a (v q ~ e q } 
x d x q + r a 


(4.4) 


The results of including this simplifying assumption are illustrated 
in figures 4.28 to 4.34. 

Neglect of the transformer voltages produces a transient re- 
sponse very similar to the original response produced with the com- 
plete alternator model. The frequency content appears to be the 
same. The maximum amplitude is increased slightly, and the maximum 
clearing time is reduced slightly. Figures 4.28 to 4.31 show the 
transient traces of the power angle for the application of a 
34 cycle three-phase short circuit fault in each model. The 
34 cycle fault is the maximum clearing time for the model neglecting 
transformer voltages . 

As a consequence of including the transformer voltages in the 
model the simulation does not indicate loss of synchronism for faults 
that cause loss of synchronism in the simpler model. The simplified 
model provides more conservative results. 

The relative values of transformer and speed voltages for the 
example studied are shown in figures 4.30 and 4.31. 
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Figure 4 , 29 ._- Comparison of rotor speed curves of baseline model with and without 

TRANSFORMER VOLTAGES IN MODEL (SAME CLEARING TIMES). 
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Wind speed 18 mph 

k d i.o 

I 1.0 

X E 0.009 PU 

Fault duration 34 cycles (0.567 sec) 
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A comparison of the swing curves for the baseline model with and 
without transformer voltages is shown in figure 4.32. These are 
swing curves for the critical clearing times for each model. The 
loss of synchronism with and without transformer voltages in the 
model are shown in figures 4.33 and 4.34 for the shortest duration 
fault causing loss of synchronism. 

4.11 Summary 

The maximum clearing time for a synchronous alternator of a hori- 
zontal axis turbine generator subjected to a three-phase, zero im- 
pedance short circuit was determined for a selected set of operating 
conditions and system parameters, namely power load, tie-line reac- 
tance, drive train stiffness, and rotary inertia. The clearing time 
determination was made by simulating a three-phase short circuit at 
the infinite bus terminals of the power system.- 

The maximum swing curve was obtained for each selected operating 
condition without including tower shadow and wind shear excitation ef- 
fects. Although these two effects are significant to the system per- 
formance of the wind turbine generator, they do not substantially af- 
fect the fault transient studied. 

The typical maximum swing curve produced by the fault applica- 
tion and removal is a positive swing of the alternator power angle of 
1 to 3 seconds duration. The two principal drive train modes of oscil- 
lation frequencies appear prominently in the transient. 

The effects upon the transient of the rotary inertia, the steady 
state power load at the time of the fault occurrence, the drive train 
stiffness, and the tie-line reactance were observed and quantified. 
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Figure 4.33, - Comparison of swing curves for baseline model with and without 

TRANSFORMER VOLTAGES SHOWING LOSS OF SYNCHRONISM. 
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Figure 4.34, - Comparison of rotor speed curves for baseline model with and without 
transformer voltages showing loss of synchronism. 
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Increased Inertia prolongs the maximum clearing time; clearing time 
is inversely ‘proportional to power load; reduced drive train stiff- 
ness tends to "decouple" the alternator and turbine; increased tie- 
line reactance reduces maximum clearing time. 



Chapter V 


CONCLUSIONS AND RECOMMENDATIONS 

The transient stability of a wind turbine generator - in par- 
ticular the response of that system to a short circuit fault - was 
investigated in this study. The scope of the problem was limited to 
a single, wind generator tied by one transmission line to an infinite 
bus. -The 100 kW DOE/NASA wind turbine generator at Plum Brook s Ohio, 
was used as the baseline system for simulation on a digital com- 
puter. 

5.1 Conclusions 
> ■ ** ^-‘ u ' - 

There are two facets to the results obtained in this research: 
one aspect is the validity of the results - how well do they reflect 
the theory; the other aspect is the magnitude of the results - what 
is the significance of the maximum clearing times calculated^ 

5.1.1 Confirmation of analysis . - The relationships observed in 
the performance of the model and listed here indicate that the re- 
sults of the simulation support the' analysis that precedes them.- 

1. Tower shadow and wind shear which produce a twice per- turbine 
revolution oscillation of the alternator power angle and thus a power 
oscillation have little or negligible effect on the transient result- 
ing from a three-phase fault. 
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2. The maximum clearing time for the wind turbine alternator sub- 
jected to a three-phase short varies directly as the square root of 
the mass moment of inertia of the system. 

3. The maximum clearing time increases as reactance between the 
alternator and the infinite bus is reduced. 

4. Maximum clearing time increases with decreasing power load. 

5. The maximum clearing time depends on the electrical load 
rather than the wind speed. When the load is doubled, the clearing 
time is halved; when the wind speed is almost doubled while the load 
remains constant, the clearing time is unchanged. 

6. As the drive train stiffness decreases so that there is less 
coupling between the turbine and the alternator, the maximum clearing 
time decreases. The relatively large inertia of the turbine becomes 
less effective in prolonging the clearing time as the coupling de- 
creases . 

The transient characteristics enumerated are in agreement with 
the power system theory discussed in Section 2 of this dissertation 
and with the phenomena observed in conventional alternating voltage 
power systems. Accordingly, the simulation model used in this research 
is valid for studying the transient phenomena of a wind turbine gen- 
erator following a three-phase short circuit fault. 

5.1.2 Significance of results . - Simulation of the wind turbine 
generator tied to an infinite bus yields maximum clearing times for 
a three-phase, zero impedance short circuit that range from 67 cycles 
(1*12 sec) to 37 cycles (0.62 sec) as the power load ranges from 
0*8 to 0*4 pu at 0,8 power factor on a "stiff" line. On a "soft" 
line the maximum clearing times range from 45 cycles (0.75 sec) to 
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21 cycles (0.35 sec) for the same load range. These times are rela- 
tively long clearing times for the average power circuit response to 
a three-phase, zero impedance fault. 

Comparison of simulation results using a model which assumes 
that the "transformer voltages" are negligible with simulation re- 
sults using a model including "transformer voltages" indicates the 
validity of assuming that the "transformer voltages" are negligible 
for the type of transient study performed for this research .- 

In order to provide a comparison between the wind turbine gener- 
ator and conventional power equipment the function plotted in fig- 
ure 5.1 is included. The figure is extracted from the General 
Electric Company publication "Electric Utility Systems and Practices" 
(28) . It is represented in that publication as characteristic of a 
"typical system." The data is presented in this discussion to il- 
lus'trate the limitation of the time scale of the switching time to 
1 second. To provide 0.-8 pu load stability the switching time for 

the ’/typical" circuit must be less than 0.1 second. This limitation 

* 

is considerably less than the 0.75 second limitation of the wind tur- 
bine generator on a "soft" line. 

The significance of the relatively long clearing time available 
with the wind turbine generator is that the system does not inpose 
a requirement for extremely fast relays and circuit breakers on its 
circuits in order to maintain transient stability. 

5.2 Recommendations for Additional Research 


The simulation and determinations of fault transients developed 


in this study lead logically to the consideration of wind turbine 





EQUIVALENT CIRCUIT ON 100,000 kVA BASE 



figure 5.1-Effec.t of switching time on maximum 
power for a three-phase fault 
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generator performance in a network. Among the problems which can 
profitably be studied with the application of the simulation and in- 
formation generated in this study are the following: 

1. The two machine problem wherein both power generators have 
the same or similar capacity. 

2 i Transient stability of a wind turbine generator tied to a bus 
by more than one transmission line. 

3. The effect of the excitation system on the stability of the 
wind power system. 

4. Transient stability for types of faults other than three-phase 
short circuits. 

5 . The effect of stability of a cluster of several wind turbine 
generators tied together electrically. 

The directions that further investigations in this field can 
take are as numerous as those which conventional power systems 
studies take. As the first studies are accomplished, how much of 
standard power system technology applies directly to wind turbine 
power systems should become progressively more evident. 



APPENDIX A 


SYMBOLS 

Alternator Constants 


L ff» L m> L JUq 
^af ’^ai&d’^aJiq 


r 


a 


x<i,xq 


r 




X 


T 1 

d 




T 


r r 
do 


T d 


rnl t 

a 

qo 


self inductance of field and damper circuits in 
' direct and quadrature axes 
Maximum mutual inductance between armature circuit 
and field or damper circuits in direct and quad- 
rature axes 

stator resistance in direct and quadrature axes 
circuits 

synchronous reactance in direct and quadrature axes 
circuits 

transient reactance in direct and quadrature axes 
circuits 

subtransient reactance in direct and quadrature axes 
circuits 

direct axis transient open circuit time constant 
direct axis subtransient open circuit time- constant 
direct axis transient short circuit time constant 
direct axis subtransient short circuit time constant 
quadrature axis subtransient open circuit time 
constant 
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v d’ v q 
v td’ v tq " 


■^d’^-q 


X d» X q 


X f 

X Ud ,X Uq 


Alternator Voltages and Flux 


bus voltage in direct and quadrature axes circuits 
stator voltages in direct and quadrature axes 
circuits 
stator voltage 

stator currents in direct and quadrature axes 
circuits 
stator current 

stator flux linkages in direct and quadrature axes 
circuits 

field flux linkage 

flux linkage in direct and quadrature axes 
e d = “ “oM a jiqi£j[,q 
. . M a£q 1 

6 d ' " “° ESq 


e ql = cooM a fif 

e q2 " “o^di-Md 


e' = 


M af 

“oEff 


M aftd 

“o L Ud X Md 


b 1> b 2 

h 

^ref 

H 


System Variables And Constants 

drive train damping factors 
height above ground 

height above ground at which wind speed = ^ Wr e f 
inertia constant 
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I 


k 


K^base- 


Ka 


K d 

Kj 

K PF 

*v 

% 




N 

P 

PF 


poles 

P„ 


P(«) 


s 


t 

T 

T 

■ L e 


mass moment of inertia 

mass moment of inertia of low speed and high speed 
components of drive train, respectively 
wind shear surface roughness coefficient 
alternator kVA rating 

spring constant equivalent of alternator synchronous 
torque 

equivalent stiffness of composite drive train spring 
gain constant of controller 

excitation gain factor for power factor control 
'partial derivative of turbine torque with respect to 
wind speed 

partial derivative of turbine torqu.e with respect to 
blade pitch 

partial derivative of turbine torque with respect to 
rotational speed 

gear ratio between high and low speed shafts 
alternator instantaneous output power 
power factor 

number of poles of alternator 
electrical power output of alternator 
damping power output of alternator 
synchronous power output of alternator 
slip of alternator; or Laplace operator 
time 

accelerating torque 
electromagnetic torque 
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■‘■in 
^sh 
T 


sync 


u 


V: 


ref 


V 


*n 


om 


V 


w 


V, 


w ref 


e 

&nom 

h 

6 


»s 


9 1» 9 2 


“o 

q A 


turbine torque 

shaft torque to alternator 

synchronous torque 

control' variable 

infinite bus voltage 

voltage regulator reference voltage 

nominal regulator voltage 

wind speed 

wind speed at reference height bygf 
state variable 

•external or tie-line reactance 
reference blade pitch 

state variables of controller and pitch servo 
alternator rotor angle 
pitch servo damping factor 

position angle of low speed and high speed compo- 
nents of drive train model 
time constant of exciter 
angular frequency of alternator 
natural frequency of pitch servo 
angular frequency of infinite bus 
angular frequency of alternator shaft 
angular frequency of turbine 



APPENDIX B 


BASELINE MODEL OF WIND TURBINE GENERATOR 

The baseline wind turbine generator system for the analytical 
model of this research is the B0E/NASA 100 kW wind turbine generator 
at the Plumbrook Station of the Lewis Research Center. This power 
generator is referred to as the Mod-0 system. The system is a pro- 
totype of the first half dozen large wind turbine systems which will 
be erected in the next several years. Although improved designs are 
expected to evolve as each system is tested, the general configura- 
tion and composition of these systems are. not expected to change 
radically. On the basis of that expectation the Mod-0 system wa$ 
esteemed to be a satisfactory model on which to base this stability 
study. A schematic diagram of the Mod-0 wind turbine generator drive 
train assembly is shown in figure B.l. 

System Description 


Configuration 

The Mod-0 machine is a horizontal axis wind turbine generator 
system mounted on top of a 30.5 meter truss tower. (General speci- 
fications of the system are listed in Table B-l.) The turbine 
(rotor) blades are downwind of the tower, an arrangement that is 
selected to provide safety from a blade striking the tower, and 
which is also more stable with respect to yawing motions. However, 
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Table B-l - General Specifications for the Mod-0 100 kW 
Wind Turbine Generator 

Horizontal axis, two-bladed down- 
wind turbine 
100 kW in 8.05 m/s wind 


Configuration 

Power Output 
Cut-in Wind Speed 
Cut-out Wind Speed 
Turbine Speed 
Yaw Rate 
Tower Height 
Tower Weight 
Blade Length 
Blade Weight 
Airfoil 
Coning Angle 
Hub 

Generator 

Weight of Components on Tower. 


4.47 m/s 
17.88 m/s 
40 r/min 
1/6 r/min 
30.48 m 
19,958 kg 
19.05 m 
907.2 kg 

NACA 23,000, ,34 deg nonlinear twist 
7 deg 

Fixed type 

Synchronous, 125 kVA, 480 v 
18,140 kg 


with this configuration the rotor blades are subject to tower shadow 
effect whereby the wind forces on a blade are abruptly and substan- 
tially reduced as the blade passes .behind the tower. 

Rotor Blades 

The rotor has two all-metal blades, each 19.0 meters long and 


each weighing 907 kg. The blades are designed to produce 133 kW of 
power at 8.05 m/s wind speed when rotating, at 40 r/min. 
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Hub 


The hub connects the blades to the low speed shaft. The hub 
also houses the mechanical gears, linkages, and actuators required 
for changing the pitch of the blades to compensate for changing wind 
and load conditions. Wind loads and centrifugal loads are absorbed 
by the hub and transmitted to the low speed shaft. The hub is a 
fixed hub; that is, the hub is bolted rigidly onto the main low 
speed shaft with the blades fixed to the hub allowing only the pitch 
change degree of freedom. 

Pitch Change Mechanism 

The pitch change mechanism consists of a hydraulic pump, a pres- 
sure control valve, an actuator, and a gear for connecting a linear 
motion of the actuator to a rotational motion (pitch) for the blades. 
The actuator is a rack and pinion that turns a master gear which in 
turn rotates the blades by means of a bevel gear mounted on the roots 
of each blade. The- hydraulic pump is mounted separately on the bed 
plate structure and the hydraulic fluid is brought into the low spe€ 
shaft- through rotating seals, thence to the pitch change components 
rotating at the other end of the shaft. 

Bed Plate 

The rotor, transmission train, alternator, and all shafts and 
bearings are mounted on a bed plate. The bed plate is supported- on 
a large gear bearing assembly which is capable of rotating (yawing) 
the entire machine on top of the tower. 
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Yaw Control 

The yaw control Is a dual control to increase structural stiff- 
ness between the tower and the rotating bed plate. Each pinion gear 
is driven by an electric motor. The yaw rate is 1/6 r/min. The con- 
trol is designed to follow slow directional changes of the wind 
rather than sudden changes. The entire structure atop the tower is 
enclosed in a fiberglass housing for protection from the environment. 

Transmission Train 

The hub transmits the high torque at low speed to the gear box 
via a low speed shaft. Out of the gear box a high speed shaft trans- 
mits the low torque at high speed, through a hydraulic slip coupling, 
to the alternator- by means of a V-belt connection. 

The gear box is a standard triple reduction design, but is used 
in a step up ratio rather than the more conventional step down ratio 
for which the gear box was originally designed. The gear box ratio 
is 45:1. 

Alternator 

The alternator is an 1800 r/min synchronous machine, three 
phase, y-connected, and self-cooled with a directly connected brush- 
less exciter and regulator.- The machine is rated at 125 kVA, 0.8 
power factor, 480 volts, 60 Hertz. 

Tower 


The tower upon which the entire wind- turbine generator is placed 
is 30.5 meters tall, constructed of steel, and is a pinned truss 



105 


design. It rests on a concrete foundation. The tower is designed to 
withstand high wind and thrust loads,- both steady and cyclic. 

Controls 

The wind turbine generator system develops a net power of 100 kW 
at wind speeds of 8.05 m/s and greater. Between 4.47 m/s and 8.05 
m/s the rotor blade pitch is fixed and the net power is a function of. 
the wind speed. At wind speeds greater than' 8.05 m/s the net power 
generated is limited to 100 kW by rotating the variable pitch- blades 
toward feather, spilling the excess energy in the wind. At wind 

4 

speeds of 17.88 m/s or greater, the turbine blades are completely 
feathered, shutting down the system. 



APPENDIX C 


DETERMINATION OF STEADY STATE INITIAL CONDITIONS 


The phasor diagram representing the initial conditions of the 
synchronous generator prior to a fault is shown in figure C.l. Given 

(1) Infinite bus voltage-, Vg 

(2) Alternator power, P, 

(3) Power factor, cos 0, 

such that- P » V t I a cos 0, determine <j>. 


AC AB + BC V B sin 4* + I a x e 


tan 9 " OA ~ QA 


V B cos <j) 


P = V fc I a COS 0 = Vgl„ COS (j) 


tan 0 = 


Vg sin <j> + x e P/V B cos $ 


Vg cos <j> 


sin 0 = JJ5 
COS 0 


Vg sin <j> cos <|> + x £ P 


> V^ cos^ 4) 


. x e P cos 0 

sin 0 cos z = sin (f> cos <j> cos 0 + ^ 


„ 2x e P cos 0 

2 sin 0 cos z <j) = 2 sin <f> cos <b cos 0 + 7 

V B 
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sin 0(2 cos 2 <{>) = (2 sin 4> cos <j>)cos 0 + 


2x e P cos 0 

V 2 

V B 


2x P P cos 0 


sin 0 (1 + cos 2(j>) = sin 2<b cos 0 + o 

V z 


B 


2x e P cos 0 

sin 0 cos 2<i> - sin 2$ cos 0 = ^2 - sin ® 

B 


2x e P cos 9 

sin(0 - 2<jj) = ^2 “ sin ® 

V B 


2x„P cos 0 


0 — 2<|> — sin - -*- ( — 6 £ 


V: 


- sin 0) 


B 


* = T 


0 + sin fsin 0 - 


2xJP cos 0 1 
© 


V: 


B 


(C.l) 


(a) 


I = 

a Vg cos <[> 


(C. 2) 


(b) 


V„ = 


Vg cos <f> 


t cos 0 


(C.3) 


= tan' 


-1 5S. 

0E 


BE = I a Xq cos <j> + I a x e cos - I a R a sin (|) 


a A e 


OE - V B + I a x e sin (j> + I a R a cos + I a x q sin <j> 
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Idxa = [j a sin(^> + fifjxd 


(d) 

V f = V B cos 5 + I aC( x d + x e ^ s: ^ n ^ + 6) + 

R a cos(<j> + 6)J 


V 

(C.5) 

(e) 

I d = I a sin (4. + «)• 

(C.6) 

(f) 

Iq = I a cos(if» + 6) 

(C . 7) 

At 

: steady state the derivatives of equations 

(3.12), (3.13), and 

(3.14) 

are zero, so that 


(8) 

E d = 0 

(C.8) 

(h) 

E q2 = 0 

(C.9) 

(i) 

E ql " V f 

(C.10) 

(j) 

A d - V - I a x« « 

(3.5) 

(k) 

A q-“ E d' " 

(3.6) 

Simultaneous solution of equations (3.5), (3.6), 

(3.9), (3.10), and 

(3.11) 

give 


a) 

E d' “ <*q - x q' , ) I q 

(C.ll) 

(m) 

' - <*d - x d> r d + V f 

(C.12) 

(n) 

V = - <X d - ’‘d'* t d + V f 

(C. 13) 
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(o) T = Ej f I, + E”I - 1. 1 (x 1 * - x*’) (3.15) 

w eddqqdq v d q ' 

Drive Train 


The initial conditions for a steady state condition of the drive 
train with a constant input torque is- given by the relationships 
shown in figure 3.4. 


^2 poles 


6 _ 

N 


(C.14) 


Turbine 


6l - 


T. - 
in 


*i®c 


K 


•D 


e 2 


(C.15) 


The correct value of 8 nom (f 3.8) for the given wind- speed 
and nominal turbine speed will produce the torque required for the 
power output desired and system losses. This value is. determined 
from the turbine torque program- (Section 3.1) and can be tabulated 
for operating conditions of interest.- For the appropriate value of 
8 n om ot h et initial conditions of the controller are zero. 



APPENDIX D 


LINEARIZED SIMULATION MODEL 

The complete set of 26 nonlinear differential equations which 
comprise the simulation model have been linearized for small per- 
turbations, and the incremental equations are listed. The block 
diagram for the system is shown in figure D.l. 

The torque developed by the turbine is- expressed as 


AT 


3T 


3T 


3T 


a AV * -t- A0, + — Ag- = K V + K,,©, + K„6. 
in 3V W w 30i 1 30 1 v w T2 1 0 K x 


R a R a x d-’ 

* d ’ f *k' - -V- + C l*d’ 5 


R x T f 
*a x q 


r: 


\ = - ~ _ t I —2, 0 I1 _ p v t t 2 

Q e q Q e d c 2 x q 6 


R a (x q x q') ( R a + x d* x q) ,, 

„ 3 — e ' + x ^ ej, - 


c 2 (x q - Xq') 6 


e ql " C 3 e q " C 4 e q' 


-q2 


°e e q 


, M+g&l ■■ . - ■■ 

« ^ n * ^3 n 


- c i c 3 (x d - X <P 6 


in 
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Figure D ♦ 1 -L i near i zed simulation model block diagram 
















113 


R s 


id “ Q e q’ + T e h' " C l 6 


K a x d 

iq = T e V " e d‘ + C 2 fi 


R a x e 


x 1 ? X 


v td = " “o~ e q' + “q - e d* + (V B cos «0 ■ C 2 x e> 6 


x fl x p R P x„ 

H e _l a e _ f i 


tq “ ~Y- e^' - (V B sin 6 0 + 0^)6 


X X 

v t ■ Z~h (- R a v tdo + x4 ? v tqo )e4’ + ~~ (x£'v tdo + RaV^)^' 

C tO^ V +*nH 


'to' 


* ^ [VB(Vtdo cos S 0 - v tqo sin 6 0 ) - x e (C 2 Vtdo + c v tqoll 6 


/ Gx” HR a \ / 

Te " [ho + Q + ~ Q/ e q' + [ho + T 1 " ”f"J e d’ " < C 1 G " G 2 H > 6 


GR a Hx^ 


el 


e q = i]Te q - 
do 




:g 2a 

T*’Q c q 

q° H 


g 3 


c 4 , , . ^ 


el + 

. e* 4- tt; — 

'J-do 

<1 

T do q T do 

K') 

t 

R a (x d - X <P 

j 

e q ■' 

T dS* 


( R i 

+ x d' x q ) ,, 


*4* + Ci 


- x r> 

rp I X 

•‘‘do 


T’ '0 
qo H 


d C 2 T" 6 
qo 


.. Bi . Kq % 

01 » - 01 "17 0 1 + if 0 2 + 


l in 


e l “ e l 
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0 2 “ 


W 


6 X *’ HR a 
[Inn + + ~^)e; 


W 


i 9 rqo - q T q ;■= q " ip \ Ido + Q 


GR a HxJ'\ 

Q J e d 


Kq Bp . Kn TJ 

+ if 01 " if 02 " if 02 * if (ClG “ c 2 h)6 


0 2 - 0 2 


j h ~ 0 2- 


e 2 = " u s e i ~ 2 ^s t 0 s 3 2 + “s 3 3 


h = - K lP ^base) 


K’ V do + R a V qo> j( C R a V do- x d’V 
Q e q + Q 


- I » 

e d 


+ [^do “• { o - v sl ” V * <°l v do - C 2'V] S 


• _ ^do x q ^-qo^a^ » t . ^do R a ~ ^gcr^d V , , t^l^do ~ Gplgp) 


i a „Q ' q 

ae^ 




where 


v e = Fieq + F 2 e‘d + F 3 V e + F 46 


v £ = v £ 


5 = 2N0r 


Q = R? + x|'x" 


V B (R a cos < 5 0 - x^' sin S Q ) 
C 1 “ Q 
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C 2 " 


V x d* cos S o + R a sin 6 o> 


_ x d ~ x d 

C 3 ~ ir' _ -sr* ' 

J x d x d 


xd - x<i 


c 4 " xi - 


x d - x d 


1 1 


K - T base (lfc- f t) 


F l-t 


k pf 

v to 


+ R 0 V„„) 
v q do a^qo 7 


. . VioK' + W 

ttWKpF i o 

ao x 


- x 


Vw , \fVtdo + K'Vtoo) 

~ + 1,1 v to Q ) 


e V 

\ to 




Kpp (R a V do x d •Vac) 


'to 


— (PF)Kpji 


(ido^a ^qo^d ^ 


*aoQ 


^o K PF _\/ X d '^tdo + ^a^tqo^ 

vi 


" x e i“2~ + 11 


'to 


VtoQ 


F~ = - — 
3 t. 


F, = 


4 t £ 


i fB (I do “» «o - I qo «*■> «o> * WiV*, - C 2''qo ) ] 
^ to 


+ (PF)K pf 


WddQ - W 


ao 


^ P o K PF 


V 


+ 1 


to 


~V V tao cos 6 o ~ V tqo sip 6 o> ~ x e (C 2 V tdo + C l V tqo 


V 


to 


® * E di - Iqo(*d’ - *q'> 
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u " T in 

z = [X dj X q > e d > e q i , e q 2 , id > iq ? vt d , v tq , v t , T J 
Matrix dimensions are x(13,l), z(ll,l)> u(l), A(13,13), B(13,l), 
C(ll,13) , K(l,24) 5 and 1(13,13)- 





MATRIX A (13x13) 
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6 


R a<*g " X <P 


R ? + X V X n 

a d g 


-C 2 <x q -xp 


C 3 (Ra + x^ T ) 

Q 


c 3 R a(*d ~ x d’> 

Q 


-Cj^Cx^ - x' d ') 


(V B cos 6 0 - C 2 x e ) 


(V B sin 6 0 + Cj^Xg) 


■V^q ( " R a V tdo + X q’ V tqo ) V^Q (X d' v tdo + VtqO 5 [?B<V> 008 5 o " V tqo ahx 5 o> “ X e< C 2 V t do + C l V t qo>] 

/ Gxq* HR a \ / ^ ^ « ♦i\ 

(V + -^ + ^) ( x + “Q T7 -<c lG -c 2 H) 


MATRIX C (11x13) 


APPENDIX E 


SYSTEM CONSTANTS 

Alternator Constants 

Base kVA = 125 kVA 

Base Voltage = 480 Volts 

Base Frequency ** 60 Hertz 

Base Ohms ~ 1..843 ohms per phase wye 


r a 

- 

0.018 

pu 


= 

2.21 

pu 


= 

0.165 

pu 


=. 

0.128 

pu 


= 

1.064 

pu 

x’ ’ 

X q 

= 

0.193 

pu 

T do 

= 

1.942 

sec 

rp t 1 

T do 

= 

0.011 

sec 

T d 

= 

0.145 

sec 

rp 1 1 

L d 

- 

0.008 

sec 

mt 1 

qo 

« 

0.062 

sec 

mt 1 

q 

= 

0.011 

sec 


Other Constants 

Bi - B2 — 1,186 • P Q N-m/.rad/sec where P^ is per unit output 
power 

h = 30.48 m 
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h ref 

= 9.14 m 

k 

* 0.22 

K I 

= 6 • 10“ 7 rad/kw-i 

k pf 

= 0.2 

N 

=• 45 

poles 

- -4 

v b 

= 1.0. pu 


= 1.0 

T o 

- 0.05 sec 

“o 

= 377 rad/sec 

“s 

= '9.-4 rad/sec • 


83 1,800 r/min 


= 4.189 rad /sec 


\ 



APPENDIX F 


SIMULATION PROGRAM LISTING 


0000100 * THREE-PHASE SHORT CIRCUIT 

0000200 

0000300 FIXED ACTUAL, COUNT , EVERY , TOTAL, COUNT R, COUNT B, COUNTH 
0000400 

0000500 FUNCTION BETICl = (1 4 . 366 , 1 . 0) , (1 8 . , O. 2439) , (30. , - 1 0 . 9878) X=O.009 

0000600 FUNCTION 3ETIC2- C 14 • 361 , 1 . 0) , ( 1 8 . , 0 . 1 847) , (30* ,-1 0 . 9733) X=0,4 

0000700 
0000800 

0000900 *XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
0001000 * X 

0001100 * INITIAL X 

0001200 INITIAL X 


0001300 *XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
0001400 

0001500 NOSORT 
0001600 

0001700 * TOTAL IS THE NUMBER OF TIMES THAT THE DYNAMIC PORTION OF THE 

0001800 * PROGRAM CYCLES. 

0001900 * COUNT IS THE MAXIMUM NUMBER OF CYCLES OF THE DYNAMIC PORTION OF THE 

0002000 * PROGRAM ALLOWED WITHOUT A TRQTUR CALCULATION. 

0002100 * KOUNTB, KOUNTR f KOUNTW ARE THE NUMBER OF CYCLES FOR WHICH THE CHANGE 

0002200 * IN BETA, RPMROT # OS WNDVEL IN ONE CYCLE OF THE DYNAMIC PORTION 

0002300 * OF THE PROGRAM IS GREATER THAN BET ALL, RP MALL, OR WNDALL, 

0002400 * RESPECTIVELY, CAUSING A NEW CALCULATION OF TRQTUR. 

0002500 


0002600 

ACTUAL 

= 

0 

0002700 

BETALL 

- 

0.005 

0002800 

COUNTR 

= 

0 

0002900 

COUNTB 

= 

0 

0003000 

COUNTW 

= 

0 

0003100 

EVERY 

= 

20 

0003200 

COUNT 


EVEFY-1 

0003300 

KOUNTR 

- 

0.0 

0003400 

KOUNTB 

= 

0.0 

0003500 

KOUNTW 

= 

0.0 

0003600 

RPMALL 

= 

0.005 

0003700 

TOTAL 

= 

0 

0003800 

TRQSET 

= 

0.0 

0003900 

TRQSUB 

= 

0.0 

0004000 

WNDALL 

= 

0.04 


0004100 
0004200 SORT 
0004300 

0004400 * ALPHA IS ANGULAR WIDTH OF TOWER SHADOW. 

0004500 ALPHA = 30. 

0004600 

0004700 ELDRAD = 62.5 

0004800 DAMPR1 = 875.*PWE?U 

0004500 DAMPS2 = DAMP31 

0005000 DAMPS = 1.0 

0005100 G3 = 6.0E-07 

0005200 INERT 1 = 95023. 

0005300 IKERT2 = 6073. 

0005400 INITAN = ICANGL* TO SAD 
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0005500 
0005600 
0005700 
0005800 
0005900 
0006000 
0006100 
0006200 
0006300 
0006400 
0006500 
0006600 
0006700 
0006800 
0006900 
0007000 
0007100 
0007200 
0007300 
0007400 
0007500 
0007600 
0007700 
0007800 
0007900 
0008000 
0008100 
0008200 
0008300 
0008400 
0008500 
0008600 
0008700 
0008800 
0008900 
0009000 * 
0009100 
0009200 
0009300 * 
0009400 
0009500 
0009600 * 
0009700 
0009800 
0009900 * 
0010000 
0010100 
0010200 * 
0010300 
0010400 
0010500 * 
0010600 
0010700 
0010800- * 


KSTIFF 

KVABAS 

NN 

OB ETA 

OH EGAS 

OMERAT 

OMEGAO 

ORPMRO 

OWNDV1 

0HNDV2 

PASSHA 

PF 

PI 

POLES 

PSI1 

PSI2 

PS 13 

PSI4 

PERIOD 

PWRNOtt 

PWRREF 

RADNOM 

BPMNOM 

RUN 

SRVLIM 

SYNSPD 

BETLIM 

TOD EG 

TOBAD 

TORPM 

TRQNOH 

TEQTOR 

WTOFPS 


3630572. 

125. 

45.0 
BETNOM 
2.*1.5*PI 
1. 

2.*PI*60. 

BPMNOH 

WHDNOM 

WNDNOH 

0.0 

0.8 

3.14159265 

4. 

0.0 

ALPHA 

180. 

ALPHA+1 80 • 

KVABAS+PHSPU 

100 . 

PHRNOM 

RPHNOM/TOBPM 

40.0 

1. 

2. ♦PI/NR 
1800. 

1.0 

180. /PI 
PI/180. 

30. /PI 

PWB N0M*NT0FPS*1 000 ./ (RADN0M*0. 8) 

(PWBPU + xA*IA*RA) *r RQNOH+DAMPR1*ICT 1D+DAHPR2*ICT2D 
0.7378 


BLOCK1 =0.0 SETS GUSAMP TO ZERO. 
BLOCK 1 = 0.0 

BLOCK 2 = 1.0 APPLIES SHORT TO LINE 
BLOCK2 =0.0 


CLTIHE IS THE TIME THAI SHORT CLEARS AND LINE IS OPENED. 
CLTXttE = 100. 


CYCLES IS DURATION OF SHORT IN 60 HERTZ CYCLES. 

CYCLES = 36. 

DURATN IS THE DURATION OF THE SHORT. 

DU RATH = CYCLES/60. 

EFLEHG = 0.0 ELIMINATES RIND SHEAR. USE 0.75 IF HIND SHEAR IS INCLUDED. 
EFLENG =0.0 

REDUCT = 0.0 ELIMINATES TOWER SHADOW. 
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.quality 


0010900 

REDUCT 

s 

o 

• 

o 

0011000 




0011100 * 

SHASRT 

IS THE TIME FOR START OF TOHEB SHADOW. 

0011200 

SHASRT 

= 

0.0 

0011300 




0011400 * 

SHOETM 

IS THE TIME AT WHICH THE SHORT IS APPLIED. 

0011500 

SHORT M 

= 

7.40 

0011600 



' 

0011700 * 

HND NOM 

IS WIND SPEED AT ZROTOR 

0011800 

WNDNOM 

- 

18. 

0011900 




0012000 * 

WNDIIH 

IS THE TIHE THAT GUST STARTS. 

0012100 

WNDTIH 

= 

1000. 

0012200 




0012300 * 



' ALTERNATOR CHARACTERISTICS 

0012400 




0012500 

RA 

= 

0.018 

0012600 

TDOP 

= 

1.94212 

0012700 

TD0PP 

ss 

0.01096 

0012800 

TQ0PP 

=J 

0.06230 

C012900 

X 

= 

0.008625 

0013000 * 

X 

ss 

0.4 

0013100 

XD 

= 

2.21+X 

0013200 

XDP 


0.165+X 

0013300 

XDPP 

= 

0.128+X 

0013400 

XQ 

= 

1.064+X 

0013500 

XQPP 

= 

0.193+X 

0013600 

XlTERM 

= 

1./X2TERH 

0013700 

X2TERM 

= 

(XD-XDPP) / (XDP-XDPP) 

0013800 




0013900 * 



VOLTAGE REGULATOR CONSTANTS 

0014000 




0014100 

KVRPF 

5= 

0.2 

0014200 

HUE 

= 

15.0 

0014300 

T1 

= 

0.05 

0014400 

VFRLIM 


3.5875 

0014500 




0014600 * 



GOST MODEL 

0014700 




0014800 

GUSPER 

=: 

1.0 

0014900 

ZREF 


30. 

0015000 

Z ROTOR 

= 

100. 

0015100 

ZZERO 

= 

0.220*3.281 

0015200 

GOSTK 

= 

3. /2. * WNDN OJI/ALOG (ZROTOR/ZZBRO) 

0015300 

GUSAMP 

= 

(GUSTK*SQBT(1.-EXP<-WNDNOH*1.467*GUSPER/. 

0015400 



(1 . 48*ZROTOF) ) ) ) *BL0CK1 

0015500 

0015600 

WNDREF 

— 

WNDNOM* (ZREF/Z ROTOR) **0.220 

v V I v* wU U 

0015700 * 



INITIAL CONDITIONS 

0015800 




0015900 * 

ICANGL 

IS AZIHUTH ANGLE OF BLADE 1 AT SYNC. 

0016000 

ICANGL 

= 

0.0 

0016100 

ICDELT 

= 

DELTA 0*TORAD/( (POLES/2. ) *NN) 

0016200 

ICBDPP 

— 

EDPP* (-TQ0 PP) 



0016300 

ICEFD 

= VF 

0016400 

ICEQP 

= EQP 

0016500 

ICEQPP 

- EQPP* (-TD0PP/X1 TERM) 

0016600 

ICLD 

= (EQPP— ID* X DPP) /OHEGAO 

0016700 

ICLQ 

= -(EDPP+IQ*XQPP) /OHEGAO 

0016800 

ICPWE 

= 0.0 

0016900 

ICS 

= 0.0 

0017000 

ICSRVR 

= 0.0 

0017100 

ICT1 

a (DAMPR1+RADNOH* (PflRPO+IA*IA*RA) *TRQNOM) /KSTIEF+ICT2 

0017200 

ICT1D 

= RADNOH 

0017300 

ICT2 

= (2. /POLES) * (DELTAR/NN) 

0017400 

ICT2D 

= ICT1D 

0017500 

0017600 PROCED PSRPU 

= PVR (WNDNOH) 

0017700 

IE (WNDNOM. LT- 18*) GO TO 1 

0017800 

PWRPU 

= 0.8 

0017900 

GO TO 

2 

0018000 1 

PWRPU 

= 0.4 

0018100 2 

CONTINUE 

0018200 BNDPRO 


0018300 

0018400 

VB 

= 1.0 

0018500 

PHI 

= ASCOS(PF) 

0018600 

CP 

= COS (PHI) 

0018700 

BETAX 

= 0. 5* (PHI+ARSIN (SIN (PHI) -2. *P*RPU*X*CP/ <VB*VB) ) ) 

0018800 

CT 

= COS (BETAX) 

0018900 

IA 

= PWRPU/ (VB*CT) 

0019000 

ARGl 

= (XQ*CT-RA*ST) *IA 

0019100 

ST 

= SIN (BETAX) 

0019200 

ARG2 

= VB+ (XQ*ST+EA*CT) *IA 

0019300 

ARG 

= ARG1/ARG2 

0019400 

DELTAS 

^ AT AN (ARG) 

0019500 

DELTA0 

= DBLTAR*TODEG 

0019600 

CD 

= COS'(DELTAR) 

0019700 

CDT 

= COS (DELTAB+BEIAX) 

0019800 

SDT 

=* SIN (DELTAR+BETAX) 

0019900 

VE 

= VB*CD+RA*IA*CDT+XD*IA*SDT 

0020000 

ID 

- IA*SDT 

0020100 

IQ 

- IA*CDT 

0020200 

EQPP 

= — (XD-XDPP) *ID+VF 

0020300 

EDPP 

= (XQ“XQPP)*IQ 

0020400 

EQP - 

= - (XD-XDP) *ID+VF 

0020500 

VD 

= VB*SIN (DELTAR) 

0020600 

VQ 

= VB*CD 

0020700 

VTD 

= VD-X+IQ 

0020800 

VTQ 

= VQ+X* ID 

0020900 

VT 

= SQRT ( VTD* VTD+VI Q*VTQ) 

0021000 

0021100 

0021200 PROCED BETNOH = TH10R2 (WNDNDM, X) 

0021300 

IF (X.EQ.0.4) GO TO 50 

0021400 

BETNOH = AEGEN (BETICl, NDNOM) 

0021500 

GO TO 

51 

0021600 50 

BETNOH = AEGEN (BETIC2, WNDN OH) 
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0021700 
0021800 
0021900 
0022000 
0022100 
0022200 
0022300 
0022400 
0022500 
0022600 
0022700 
0022800 
0022900 
0023000 
0023100 
0023200 
0023300 
0023400 
0023500 
0023600 
0023700 
0023800 
0023900 
0024000 
0024100 
0024200 
0024300 
0024400 
0024500 
0024600 
0024700 
0024800 
0024900 
0025000 
0025100 
0025200 
0025300 
00254 00 
0025500 
0025600 
0025700 
0025800 
0025900 
0026000 
0026100 
0026200 
0026300 
0026400 
0026500 
0026600 
00267 00 
0026800 
0026900 
0027000 


ORIGINAE PA.GS1S 
OF POOR QUALHS 


51 CONTINUE 

ENDPRO 


VFO 

= VF 

VNOM 

= VT 

VRZ 

= 0.0 

VTDEL 

= 0.0 


* DEBUG 

NO SORT 

1000 CONTINUE 

CALL DEBUG (1 #0.0) 

♦xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

* 

* DYNAMIC 

DYNAMIC 

♦xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


* 


.DRIVE TRAIN. 


TH1 

OMEGA 1 

HPMROT 

TRQTPU 

TH1DD 

TH2 

OMEG&2 

SPEED 

ALTRPM 

TH2DD 

RADDPS 

RPMERR 

OMEGAS 

DEL TAR 


INTGRL (ICT1,OMEGA1 ) 

INTGRL (ICT1D/IH1DD) 

OMEGA1*TORPH 
TRQIUR/TRQ NOM 

(TRQTUB-KSTIFF* (TH1-TH2) —DAMPR1 IOMEGA 1 ) /INERT1 
INTGRL (ICT2,OMEGA2) 

INTGRL (ICT2D,TH2DD) 

ALTRPM/SYNSPD 

OMEGA2*TORPH*NN 

(KSTIFF*(TH1-TH2) - DAHPR2*OMEGA2-TRQLOD) /INERT2 

0HEGA2— RADNOH 

2ADDPS*T0RPM 

OHEGAO+ {POLES/2- ) *RADDPS *NN 
(POLES/2- ) *NN*INTGRL (ICDELT , RADDPS) 


DELTA D = DELTAR*TODBG 
ABDELT = ABS(DELTAD) 


♦ 


CONTROLLER AND SERVO. 


PHRERR 

SERVIN 

SBVERR 

LI MIN 

AN GIN 

LIMOUT 

AN GOUT 

SERVOT 

BETA1 

BETERR 

BETODT 


PHRBEF-PWRLO0 

G3* INTGRL (ICPHR,PffRERR) 

OMEGAS+OHEGAS* 5 (SERVIN- SERVOT) -2. *DA»PS*OMEGAS*LIKIN 
INTGRL (IC5RVR, S RVERR) 

LIMIN*T0DEG 

LIMIT (-SRVLlMf SRVLIM, LIMIN) 

LIMOUT*TODEG 
INTGRL (ICS, LIMOUT) 

SERVOT 

BETA1*T0DEG 

BETERR+BETNOM 
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BETA LIMIT 


0027100 * 

0027200 

0027300 PSOCED BETA = TH£T(BETOUT,BETLIM) 

0027400 BETA = BETOUT 

0027500 IF (BETA. LE. BETLIM) GO TO 40 

0027600 BETA = BETLIM 

0027700 40 CONTINUE 

0027800 ENDPRO 

0027900 

0028000 * TURBINE TORQUE DEVELOPED 

0028100 

0028200 PROCED TRQTUR, TRQTU1 ,TRQTU2, ACTUAL, COUNT, TOTAL, KOUNT B, KOUNTR, KOUNTW, .. - 
0028300 COUNTB, COUNTR, COUNTW = TRQ (WNDVE1 , WNDVE2 ,RPHROT, BETA, EVERY,, . . 

0028400 OBETA,ORPHRO,OWNDV1 , OWNDV2 , BET ALL, RPMALL, WNDALL, TRQSUB) 

0028500 TOTAL = TOTAL+1 

0028600 COUNT = COUNT+1 

0028700 TRQSET = 0.0 

0028800 IF (COUNT. EQ. EVERY) GO TO 14 

0028900 BETABS = ABS (BETA-OBETA) 

0029000 RPMABS = ABS (RPMROT-ORPMRO) 

0029100 WNDAB1 = ABS { HNDVE1-OWNDV1 ) 

0029200 WNDAB2 = ABS (WNDVE2-OWNDV2) 

0029300 IF (BETABS »LE« BETALL) GO TO 11 

0029400 COONTB - COUNTB+1 

0029500 KOUNTB = COUNTB 

0029600 TRQSET =1.0 

0029700 11 IF (RPKABS.LE. RPMALL) GO TO 12^ 

0029800 COUNTR = COUNTTl+1 

0029900 KOUNTB = COUNTR 

0030000 ' TRQSET =1.0 

0030100 12 IF (WNDAB1 • LE. WNDALL) GO TO 13 

0030200 IF (WNDAB2.LE. WNDALL) GO TO 13 

0030300 COUNTW = CODNTW+1 

0030400 KOUNTW = COUNTW 

0030500 TRQSET =1.0 

0030600 13 IF (TRQSET. EQ. 1 . 0) GO TO 14 

0030700 TRQTUR = TSQSUB 

0030800 GO TO 15 

0030900 14 TRQ1,THH1,BET1 = TORKZ { WNDVE1 , RPMROT,BETA) 

0031000 TSQTU1 = 0.5*TRQl 

0031100 TRQ2,THR2,BET2 = TORKZ (WNDVE2 , RPMROT,BETA) 

0031200 TRQTU2 = 0.5*TBQ2 

0031300 TRQTUR = TRQTU1+TRQTD2 

0031400 ACTUAL = ACTUA1+1 

0031500 TRQSUB = TRQTUE 

0031600 COUNT = 0 

0031700 15 OBETA = BETA 

0031800 ORPMRO = RPMROT 

0031900 OHNDV1 = HNDVE1 

0032000 OWNDV2 = WNDVE2 

0032100 16 CONTINUE 

0032200 ENDPRO 
0032300 

0032400 * _ TOWER SHADOW . 
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0032500 

0032600 

0032700 

0032800 

0032900 

0033000 

0033100 

0033200 

0033300 

0033400 

0033500 

0033600 

0033700 

0033800 

0033900 

0034000 

0034100 

0034200 

0034300 

0034400 

0034500 

0034600 

0034700 

0034800 

0034900 

0035000 

0035100 

0035200 

0035300 

0035400 

0035500 

0035600 

0035700 

0035800 

0035900 

0036000 

0036100 

0036200 

0036300 

0036400 

0036500 

0036600 

0036700 

0036800 

0036900 

0037000 

0037100 

0037200 

0037300 

0037400 

0037500 

0037600 

0037700 

0037800 


ORIGINAL PAGE IS 
©S' POOR QUALITY! 


PROCED REHAIN,SHAD01,SHAD02, GAMDEG. = SHADO (RPMROT,PASSHA, SHASRT) 
SHAD01 =0,0 
SHAD02 = OlO 

PSIBAD « INITAN + O .O/TOHPM) *INTGBL (0 . 0,RPMR0I) 

PSIDEG = PSIRAD*TODEG 
REMAIN = AHOD (PSIDEG, 360.) 

ANGLE = REMAIN 

IF ((REMAIN.GE.PSI1) .AND. (REMAIN. IE. PSI2)) ANGLE = PSI1 
IF ( (REMAIN. GE. PSI3) * AND. (REMAIN. IT. 210. ) ) ANGLE = PSI-3 
GAMDEG = (180./PSI2)* (REMAIN-ANGLE) 

IF {PASSHA.EQ. 1.0) GO TO 22 

IF ( (ABS (GAMDEG) ,LE. 0.001) .AND. (TIME. GE. SHASRT) -AND. ... 

{REDUCT. NE. 0.0) ) GO TO 22 
GO TO 21 

22 PASSHA = 1.0 

GAMRAD = GAMDEG*! ORAD 

IF (ANGLE * EQ- PSI3) GO TO 20 

SHAD01 = RED0CT*SIN (GAMRAD) 

GO TO 21 

20 SHAD02 = REDOCT*SIN (GAMRAD) 

21 CONTINUE 
ENDPRO 


SHADOW = SHAD01 +SHAD02 

* VOLTAGE REGULATOR 

VF A MINI (VFO + VFDEL, VFRLIM) VI 

VTDEL = { VNOM- VT) - KVEP F* (PF*IA- ( VTD*ID+ VTQ*IQ ) /VT) 

ZINTER = REALPL (0- 0,,T1, VTDEL) 

VRZ = MUE*ZINTER 
VFDEL = INTGRL (0. 0 ,VRZ) 

* ALTERNATOR EQUATIONS 

PROCED OMERAT = OR ATIO (BLOCK2 , OMBGAO, SHORTM, OMEGAS) BL0CK2 

IF (BL0CK2. EQ. 1 ■ 0} GO TO 70 BLOCK2 

71 OMERAT =1.0 
GO TO 72 

7 0 IF (TIME. LT. SHORT M) GO TO 71 

OMERAT = OMEGAF/OHEGAO 

72 CONTINUE 

ENDPRO 

PROCED VD,VQ = VOLTS (DELTAR, OMERAT, SHORTM ,DURATN,CLTIME , EDPP, EQPP) 

IP ( (TIME. GE.CLTXME) .AND. (TIME. LT. (SHORTM+DURATN) ) ) GOTO 30 

SWITCH = 1.-BLOCK2* (STEP (SHORTM) -STEP (SHORTM+DURATN) ) BLOCK2 

VD = SIN(DELTAR) *SVITCH 

VQ = COS (DELTAS) ♦SWITCH 

GO TO 31 

3 0 VD = ED PP* OMERAT 

VQ =» EQPP+OMERAT 

31 CONTINUE 
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0037900 ENDPRO 
0038000 

0038100 VTD = VD-X*IQ 

0038200 VTQ = VQ+X+ID 

0038300 VT = SQRT ( VTD* VTD+VTQ* VTQ) 

0C384 00 

0038500 PROCED PFLOAD = LO ADPF (VT , VTD, VTQ , IA, ID ,IQ) 

0038600 IF ((VT*IA) * EQ. 0.0) GO TO BO 

0038700 P FIORD = (VTD* ID+VTQ*IQ) /( VT*IA) 

0038800 GO TO 81 

0038900 80 P FLO AD = 0.0 

0039000 81 CONTINUE 

0039100 ENDPRO 
0039200 

0039300 TERMLD = ( VD+RA*ID- (ED+XQ*IQ) *OMERAT) 

0039400 TERMLQ = (VQ+RA*IQ* (XD*ID-EQ1-EQ2) *OHEKAT) 

0039500 LAMDAD = OMEGA 0*INTGRL (ICLD , TERMED) 

0039600 LAMDAQ = OMEGAO*INTGRL(ICLQ, TERMLQ) 

0039700 LMDDOT = DERIV (0 .0, LAMDAD) 

0039800 LMQDOT = DERIV (0 .0, LAMDAQ) 

0039900 RATIO D = LMDDOT/ (OMEGAO*! AM DAD) 

0040000 RATIO Q = LMQDOT/ (OMEGA0*LAMDAQ) 

0040100 EDPP = (-1 ./TQOPP) +INTGRL (ICEDPP, ED) 

0040200 EQPP = (-X1TERM/TD0PP) *INTGRL(ICEQPP,EQ2) 

0040300 TER Ml = (VF-BQ1) /TDOP 

0040400 EQP * INTGRL (ICBQP,TERM1 ) 

0040500 ID = (EQPP- LAMDAD) /XDPP 

0040600 EQ1 = EQP*X2TEBH-EQPP* (XD-XDP) /(XDP-XDPP) 

0040700 EQ2 = ( ((XDP-XDPP) *LAHDAD+XDPP*EQF-XDP*EQPP) *X2 TERM/ (-XDPP) ) 

0040 BOO IQ = (EDPP+LAMD AQ) / (-XQPP) 

0040900 ED = < (LAMDAQ* (XQ-XQPP) +EDPP*XQ)/XQPP) 

0041000 IA = SQRT (ID*ID+LQ*IQ) 

0041100 TALTPU = EQPP*IQ+E DPP* ID-ID* IQ* (XDPP-XQPP) 

0041200 TRQALT = TALTPU*TRQNOM/NN 

0041300 TRQ10D = TBQALT*NN 

0041400 PWRPU = VD*ID+VQ*IQ 

0041500 PWBLOD = KVABAS*PHRPU 

0041600 TORQIN as TBQTUR-DflflPRI+OHBGA 1-DAMPR2*OHEGA2 

0041700 PWRIN = T0RQIN*0«EGA2/ (1000. *WTOFPS) 

0041800 

0041900 * TORQOE LOSS 

0042000 

0042100 TRQLOS = DAHPRl* (OMEGA1 +OHEGA2) +INERT1*TH1 DD+INBRT 2*TH2DD. . . 

0042200 +INERT2+TH2DD 

0042300 

0042400 * WIND SPEED _ 

0042500 

0042600 ZBLAD1 « ZROTOR+BLDRAD*EFLENG*COS ( (REMAIN-ALPH A/2 . ) *TORAD) 

0042700 ZBLAD2 = ZROTOR-BLDRAD*EFLENG*COS ( (REMAIN- ALPHA/2. ) *T0RAD) 

0042800 WNDH01 = WNDREF* (ZBLAD1 /ZREP) *+0.220 

0042900 WNDN02 « WNDREF*, (ZBLAD2/ZREF) **0.220 

0043000 WNDVE1 = (WNDN01 +GUST) * (l.-SHADOl) 

0043100 WNDVE2 = (HNDN02+GUST) * ( 1 . -SHAD02) 

0043200 WNDAVE = 0 .5* (WNDVE1 +WNDVE2) 
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0043300 GOST = GUSAMP* (1 . -COS (2- *PI* (TIME-WNDTIM) /GOSPER) ) . . . 

0043400 * (STEP (KN DTI ft) - STEP (HNDTIM+GUSPER) ) 

0043500 PERCNT “ 1 00 . *GUST/WNDNOM 

0043600 

0043700 PROCED ILIHE = LIKEI (ID,IQ, OMEGAS, TIME) 

0043800 IF(TIME.LT. (SHORTM-O.5) ) GO TO 90 

0043900 ILIHE = ID*COS (OMEGAR*! IHE) ~IQ*S IH (OMEGAR*TIHE) 

0044000 90 CONTINUE 

0044100 ENDPRO 
0044200 

0044300 *XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
0044400 * 

0044500 * ' TERMINAL 

0044600 TERMINAL 


0044700 

0044800 

0044900 

0045000 

0045100 

0045200 

0045300 

0045400 

0045500 

0045600 

0045700 

0045800 

0045900 

0046000 

0046100 

0046200 


♦xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


200 


100 


WRITE (3,200) DELTAD, TIME, CYCLES 
WRITE (6, 200) DELTAD, TIME, CYCLES 

FORMAT C///5X , 'DELI AD- *,F7.2,5X, • TIMB= F7. 3, 5X , 1 CYCLES- F4- 0///) 
WRITE (06,100) EVERY, TOTAL , ACTUAL, G3, PWSPU, KSTIPF ,WNDNOM, RUN, BA, • . . 

DB1TA0 ,X, BET NOH,BPLENG, BLOCK1, BLOCK2, REDUCT, DAHPR1 , CYCLES 
FORMAT {1H0 ,7 HEVE1Y = ,I3,8X, 7HTOTAL =, 16 , 8X, 7HACTUAL=, 16 , 

$ 8X,7HG3 =,F10.7,2X,7HPOWER = ,F4. 1 ,7X,7HKSTIFF=,F10. 0//'- 

$ 1X,7HWNDNOH=,F7-3,4X,7HRUN =,F6. 0, 8X, 7HRA =,F5.3,9X, 

$ 7HDELTA0=,F8.3,4X,7HX F6. 3 ,5X,7HBETN0M=,F1 1 . 6// 

$ 1X,7HEFLENG=,F5-2,6X,7HBL0CK1=,F4. 1 , 10X , 7HBLOCK2=, F4 . 1 , 

$ 10r # 7HREDUCT=,F4.2,8X,7HDAMPR1=,P5.0,6X,7HCYCLES=,F4.0) 


METHOD RKSFX 


0046300 RANGE DELTAD, SPEED, K OUNTB, KOUH TR, KOUNTW 
0046400 

0046500 PRTPLT DELTAD (SPEED, TBQTUR, PWRPU ) 

0046600 

0046700 PREPARE DELTAD, SPEED, PWRPU, TRQTPU, SHADOW, WND7E1 
0046800 

0046900 FINISH ABDELT « 200. 

0047000 

0047100 TIMER FINTIM = 2.0, OUTDEL = 0.002, DELT = 0.002 
0047200 
0047300 END 
0047400 STOP 
0047500 ENDJOB 
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